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FOREWORD 


The  work  described  in  this  Annual  Technical  Report  was  performed  at  the 
United  Technologies  Research  Center  (UTRC)  under  Contract  N00014-79-C-0344 , 
Modification  P00002,  entitled  "Parametric  Analysis  and  Optimization  of  Two- 
Phase  Nozzles",  for  the  Office  of  Naval  Research  (ONR) .  This  report  summarizes 
results  obtained  for  the  Phase  II  -  (second  year)  study  indicated  above  which 
is  an  integral  part  of  a  three-year  study  program  on  Two-Phase  Gas-Liquid 
Nozzle  Flow  Theory.  Dr.  Simion  C.  Kuo  is  the  Principal  Investigator  for  this 
contract,  and  Dr.  C.  W.  Deane  is  the  major  contributor  for  this  phase  of  the 
study. 

The  contract  program  was  initiated  with  ONR  on  May  15,  1980,  and  the  ONR 
Program  Manager  is  Mr.  M.  Keith  El lingsworth,  Power  Program,  ONR,  Arlington, 
Virginia.  Valuable  guidance  and  comments  received  from  Mr.  Ellingsworth  are 
gratefully  acknowledged. 
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Two-Phase  Nozzle  Theory  and  Parametric  Analysis 
Phase  II  -  Parametric  Analysis  and  Optimization 


SUMMARY 


The  objectives  of  the  work  presented  in  this  report  were  to  analyze  the 
parametric  performance  and  to  identify  the  optimum  design  criteria  for  two-phase 
nozzles.  A  working  theory  of  two-component  two-phase  (gas-liquid)  nozzles  developed 
earlier  (as  part  of  the  Phase  1  work  on  this  program)  was  used  to  conduct  these 
parametric  analyses  over  a  wide  range  of  flow  conditions  and  fluid  properties; 
this  working  theory  includes  droplet  breakup  and  the  interphase  transport  of 
momentum  and  thermal  energy  but  excludes  mass  transfer  and  frictional  losses.  The 
sensitivity  of  the  predicted  nozzle  performance  to  the  assumptions  concerning  the 
droplet  characteristics  was  examined,  and  criteria  were  developed  for  optimum 
nozzle  designs  based  on  the  parametric  analyses.  A  simplified  procedure  was 
developed  for  the  purpose  of  estimating  the  approximate  dimensions  of  two-phase 
nozzles.  The  performance  criteria  of  two-phase  nozzles  were  then  assessed  for 
the  nozzle  as  an  integral  part  of  an  energy  conversion  system. 

The  sensitivity  of  the  predicted  nozzle  performance  to  the  droplet  character¬ 
istics  was  examined  by  varying  the  correlations  assumed  in  the  formulation  of  the 
working  theory  and  then  estimating  the  resulting  changes  in  nozzle  performance. 

Then,  using  the  Model  with  Droplet  Heat  Transfer,  parametric  performance  analyses 
of  two-phase  nozzles  were  conducted  as  a  function  of  the  governing  flow  parameters, 
including  loading  ratio,  inlet  and  exit  pressures,  inlet  temperatures  of  the  gas 
and  of  the  liquid,  thermophys ica 1  properties  of  the  liquid  and  of  the  gas,  inlet 
velocities,  droplet  size,  and  nozzle  length.  Efficiency  maps  of  a  two-phase 
nozzle  were  then  constructed  based  on  the  results  of  these  analyses,  and  the  power 
split  between  the  liquid  and  the  gas  streams  was  estimated  under  different  operating 
condi t ions . 

Based  on  the  parametric  analyses,  criteria  for  optimum  nozzle  designs  were 
developed,  and  desirable  liquid/gas  combinations  were  investigated  in  terms  of 
three  criteria:  maximum  operating  temperatures  of  the  liquid,  degradation  charac¬ 
teristics  of  the  liquid,  and  desirable  thermophysical  properties.  A  simplified 
nozzle  sizing  procedure,  based  on  isothermal  theory,  was  devised  for  the  quick 
calculation  of  the  approximate  dimensions  of  the  nozzle.  Finally,  the  work 
extraction  potential  of  two-phase  nozzles  was  examined  by  identifying  nozzle 
performance  criteria  and  by  assessing  the  attributes  of  a  nozzle  suitable  for 
incorporation  in  an  energy  conversion  system  in  terms  of  the  off-design  performance 
characterist ics  of  the  nozzle. 
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This  study  program  was  conducted  by  the  Energy  Conversion  Systems  Analysis 
Group  at  UTRC  under  Contract  N00014-79-C-0344 ,  Modification  P00002,  from  the 
Office  of  Naval  Research,  Power  Program,  Arlington,  Virginia. 
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INTRODUCTION 


Power  conversion  systems  involving  two-pha6e  working  fluids  are  of  potential 
interest  in  a  number  of  applications  because  of  several  inherent  advantages.  For 
example,  two-phase  work  extraction  can  be  directly  applied  to  harness  geotherma’ 
energy,  where  the  working  fluid  is  initially  available  as  a  two-phase  mixture  of 
steam  and  brine.  Another  potentially  attractive  application  is  in  the  marine 
propulsion  field  where  the  low  exit  velocities  of  a  high-density  two-phase  mixture 
from  a  two-phase  nozzle  (coupled  with  a  two-phase  turbine  and  other  components  to 
form  a  work  extraction  system)  can  advantageously  drive  the  propulsion  turbines 
with  shaft  speeds  significantly  lower  than  those  in  more  conventional  turbomachinery 
(steam  or  gas  turbine)  power  systems.  Lower  shaft  speeds,  in  turn,  would  allow  the 
typically  massive  gearbox  (that  couples  the  turbine  shaft  to  the  thrustor  shaft) 
to  be  much  smaller  and  lighter  (because  less  speed  reduction  would  be  required), 
thereby  allowing  increased  payload  or  range.  Another  potential  advantage  in 
marine  propulsion  applications  is  the  possible  elimination  of  the  boiler.  This 
might  be  possible  if  the  liquid  phase  (for  example,  a  heat  transfer  oil  with  low 
vapor  pressure)  were  heated,  under  pressure  without  boiling,  and  then  used  to 
flash  the  condensate  of  the  second  component  (for  example,  water)  into  the  gas 
phase  (steam  in  this  example). 

In  general,  a  two-phase  (gas-liquid)  mixture  can  contain  the  same  amount  of 
enthalpy  at  a  much  lower  stagnation  temperature  than  a  single-phase  gas  such  as 
steam.  Because  of  the  relatively  high  density  of  the  liquid,  the  two-phase 
mixture  can  also  carry  the  same  amount  of  kinetic  energy  in  the  nozzle  exit  stream 
at  a  much  lower  velocity  than  the  single-phase  gas.  For  example,  the  differences 
in  velocity  and  temperature  levels  between  the  single-phase  (l-$)  and  the  two-phase 
(2-$)  streams  are  indicated  in  Fig.  II. I  where  both  streams  have  the  same  gas 
mass  flow  rate  (of  steam),  but  the  two-phase  stream  contains  liquid.  For  equal 
amounts  of  exit  kinetic  energy  and  at  loading  rates  greater  than  ten,  the  exit 
velocities  of  the  two-phase  mixture  are  generally  less  than  one-third  of  the  exit 
velocities  of  the  single-phase  system.  These  lower  nozzle  exit  velocities  then 
allow  lower  turbine  shaft  speeds.  Further,  to  achieve  the  same  exit  gas  temperature, 
the  nozzle  inlet  temperatures  of  the  single-phase  system  with  steam  are  generally 
a  minimum  of  350  R  higher  than  the  inlet  temperatures  of  a  two-phase  system  with 
steam  and  Dowtherm  A  (DTA) ,  for  example.  Hence  the  two-phase  system  also  has  the 
potential  to  operate  with  lower  heat  source  temperatures. 

In  addition  to  the  two  applications  cited  above  for  two-phase  work  extraction 
systems,  many  other  applications  involve  flows  of  a  substance  that  is  a  suspension 
of  liquid  droplets  in  a  gas.  The  particular  area  of  interest  in  this  investigation 
is  a  two-phase  nozzle  where  a  gas  and  a  liquid  at  conditions  of  high  pressure  and 
low  velocity  are  mixed  at  the  nozzle  inlet  and  then  expanded  through  the  nozzle  to 
conditions  of  low  pressure  and  high  velocity.  Previous  investigators  have  proposed 
physical  models  for  two-phase  flows  in  nozzles,  but  these  models  generally  display  one 
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of  two  drawbacks  in  terms  of  conducting  parametric  analyses.  Either:  only  some 
of  the  significant  effects  of  momentum  and  heat  transfer  were  considered  (e.g., 

Refs.  1  to  7);  or  these  effects  were  considered  but  the  model  (e.g.,  Refs.  8  and 
9)  was  formulated  essentially  for  the  purpose  of  detailed  theoretical  computation  and 
not  for  the  purpose  of  conducting  parametric  analyses  and  optimization  where  many 
cases  must  be  examined  quickly  yet  with  reasonable  accuracy.  Hence,  to  avoid  these 
drawbacks,  a  more  practical  working  theory  was  formulated  during  Phase  I  of  this  study 
program  (Ref.  10).  This  working  theory,  developed  and  computerized  as  a  Model  with 
Droplet  Heat  Transfer,  is  suitable  for  parametric  analyses  and  optimization  and  is 
formulated  for  a  wide  range  of  operating  conditions  because  it  includes  the  significant 
effects  of  velocity  slip,  droplet  breakup,  and  droplet  heat  transfer. 

The  Model  with  Droplet  Heat  Transfer  assumes  a  two-component,  two-phase 
mixture,  where  the  liquid  phase  is  in  the  form  of  droplets.  One-dimensional  flow, 
with  variables  changing  only  in  the  axial  flow  direction  of  the  nozzle,  was 
selected  for  study  so  that  both  the  model  and  the  results  will  be  tractable. 

Numerical  marching-type  solutions  with  the  nozzle  divided  into  numerous  small 
axial  segments  of  equal  length  were  chosen  because  closed-form  analytical  solutions 
are  possible  only  for  certain  simplified  situations.  The  variations  of  the 
thermophysical  properties  of  the  gas  and  the  liquid  with  temperature  and/or 
pressure  along  the  nozzle  are  included  in  the  marching-type  solution  by  recalculat¬ 
ing  the  value  of  the  property  for  the  local  temperature  or  pressure  at  each 
segment . 

The  Model  with  Droplet  Heat  Transfer  uses  the  following  six  basic  equations 
to  describe  the  adiabatic  one-dimensional  two-phase  nozzle  flow: 

Cont inui ty 

Momentum 

Droplet  Drag 

Droplet  Breakup 

Energy 

Heat  Transfer 

These  equations  describe  the  flow  conditions  in  and  across  each  segment  of  finite 
length,  and  each  equation  has  been  transformed  into  a  form  suitable  for  computerized 
numerical  solutions,  as  described  in  detail  in  Ref.  10. 

The  objectives  of  the  Phase  II  -  Two-Phase  Nozzle  study  program  were  to  conduct 
parametric  analyses  using  the  Model  with  Droplet  Heat  Transfer  that  was  developed 
during  Phase  I,  and  to  identify  optimum  design  criteria  for  two-phase  nozzles.  The 
results  of  this  Phase  II  study  program  are  presented  in  this  report. The  sensitivity 
of  the  two-phase  nozzle  performance  (Model  with  Droplet  Heat  Transfer)  to  the 
assumptions  concerning  the  droplet  characteristics  was  analyzed  in  Task  II. 1  by  varying 
the  correlations  assumed  in  the  formulation  of  the  theory.  Then  this  theory  was 
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used  to  conduct  parametric  performance  analyses  in  Task  11.2  by  examining  the 
effect  of  the  governing  flow  parameters  on  nozzle  performance  predictions  and 
efficiency,  including  the  effects  of  loading  ratio,  nozzle  pressure  ratio,  nozzle 
length,  liquid  surface  tension  and  the  related  droplet  size,  as  well  as  the 
effects  of  other  flow  parameters.  Based  on  these  parametric  analyses,  criteria 
were  developed  in  Task  II. 3  for  optimum  nozzle  designs;  desirable  liquid/gas 
combinations  were  investigated,  and  a  simplified  procedure  for  two-phase  nozzle 
sizing  was  developed.  Finally,  two-phase  nozzle  performance  criteria  applicable 
to  work  extraction  were  identified  in  Task  II. 4,  and  a  simplified  procedure  was 
developed  for  the  calculation  of  off-design  nozzle  performance  characteristics. 


R81-955229-4 


INTRODUCTION  REFERENCES 


1.  Tangren,  R.  F. ,  C.  H.  Dodge  and  H.  S.  Siefert:  Compressibility  Effects 
in  Two-Phase  Flow.  J.  Appl .  Phys.,  Vol.  20,  No.  7,  pp.  637-645,  1949. 

2.  Kliegel,  J.  R. :  One  Dimensional  Flow  of  a  Gas-Particle  System.  Report 
No.  TR- 59-0000-00746 ,  Space  Technology  Laboratories,  Los  Angeles,  1959. 

3.  Rudinger,  G.:  "Relaxation  in  Gas-Particle  Flow,"  Chapter  in  Nonequilibrium 
Flows  —  Part  I,  P.  P.  Wegener,  ed.,  pp.  119-161.  Marcel  Dekker,  Inc., 

1969. 

4.  Netzer,  D.  W. :  Calculations  of  Flow  Characteristics  for  Two-Phase  Flow 

in  Annular  Converging-Diverging  Nozzles.  Report  No.  TM-62-3,  Jet  Propulsion 
Center,  School  of  Mechanical  Engineering,  Purdue  University,  1962. 

5.  Hultberg,  J.  A.  and  S.  L.  Soo:  Two-Phase  Flow  Through  a  Nozzle.  Astronautic 
Acta,  Vol.  11,  No.  3,  pp.  207-216,  1965. 

6.  Crowe,  C.  T.,  M.  P.  Sharma  and  D.  E.  Stock:  The  Partic le-Source-In  Cell 
(PSI-Cell)  Model  for  Gas  Droplet  Flows.  ASME  J.  of  Fluids  Eng’g.,  pp. 
325-332,  1977. 

7.  Comfort,  W.  J.,  T.  W.  Alger,  W.  H.  Giedt  and  C.  T.  Crowe:  Calculation  of 

Two-Phase  Dispersed  Droplet-in-Vapor  Flows  Including  Normal  Shock  Waves. 

ASME  J.  of  Fluids  Eng'g.,  pp.  355-362,  1978. 

8.  Elliott,  D.  G.:  "Theoretical  and  Experimental  Investigation  of  a  Gas-Driven 
Jet  Pump  for  Rocket  Engines",  in  Liquid  Rockets  and  Propellants,  Progress  in 
Astronautics  and  Rocketry  -  Vol.  2,  pp.  497-541.  Academic  Press,  1960. 

9.  Elliott,  D.  G.  and  E.  Weinberg:  Acceleration  of  Liquids  in  Two-Phase  Nozzles 
JPL  Report  32-987,  1968. 

10.  Deane,  C.  W.  and  S.  C.  Kuo:  Two-Phase  Nozzle  Theory  and  Parametric  Analysis, 
UTRC  Report  R80-954624-2 ,  1980. 


6 


R81-955229-4 


CONCLUSIONS 


1. 
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n  percent)  of  a  two-phase  two-component  nozzle 
d  60' s  and  the  high  80's  depending  on  the  loading 
quid-to-gas  mass  flowrates),  the  nozzle  pressure 
and  to  a  lesser  degree  the  inlet  conditions  and 
ds.  In  contrast,  a  properly  designed  single- 
straight  axis  can  achieve  efficiencies  ranging 


2.  Generally  speaking,  the  higher  end  of  the  efficiency  range  indicated  above 
can  be  achieved  in  a  long  nozzle  (say  greater  than  10  inches)  operating  at 
a  relatively  low  pressure  ratio  (say  under  50)  while  the  lower  efficiencv 
range  is  for  shorter  nozzles  particularly  with  a  very  high  pressure  ratio 
(up  to  several  hundreds).  Using  a  mixture  of  steam  and  Dowtherm  A(DTA), 
a  six  inch  nozzle  operating  at  a  loading  ratio  of  40  and  a  pressure  ratio 
of  40  would  have  an  efficiency  between  70  and  75  percent. 


3.  At  a  given  pressure  ratio,  the  nozzle  efficiency  goes  through  a  minimum  as 

the  loading  ratio  is  increased  —  at  first  decreasing  because  of  the  presence 
of  more  liquid,  reaching  a  minimum  value  (at  a  loading  ratio  in  the  vicinity 
of  10),  before  increasing  gradually  because  the  slip  velocity  between  the 
two  phases  is  decreasing  as  the  loading  ratio  is  increased  to  approximately 
60.  At  higher  loading  ratios,  the  flow  regime  may  no  longer  be  droplet 
flow,  as  assumed  in  the  working  theory. 


4.  The  exit  specific  power  (flux  of  kinetic  energy  at  the  nozzle  exit)  decreases 
when  either  the  loading  ratio  or  the  pressure  ratio  or  both  is  increased. 

As  the  loading  ratio  is  increased  (more  liquid),  the  required  nozzle  exit 
area  also  increases  because  the  gas  exit  velocity  is  decreasing  and  more 
flow  area  is  needed  to  pass  the  lower-velocity  mixture. 

5.  For  a  mixture  of  steam  and  Dowtherm  A(DTA)  flowing  at  a  loading  ratio  of  40, 
the  power  split  (between  the  liquid-phase  kinetic  energy  and  the  gas-phase 
kinetic  energv  at  the  nozzle  exit)  is  24  for  a  nozzle  operating  at  an  inlet 
pressure  of  1^0n  nsia  and  a  pressure  ratio  of  10.  With  the  pressure  ratio 
increased  to  40  for  these  inlet  conditions,  the  power  split  is  20.  Further, 
for  loading  ratios  greater  than  five  for  steam/DTA  mixtures  at  those  operat¬ 
ing  conditions,  over  75  percent  of  the  mixture  kinetic  energy  is  carried  by 
the  liquid  phase. 


6.  The  sensitivity  of  the  predicted  nozzle  performance  to  assumptions  of  droplet 
breakup  was  found  to  be  significant  (as  much  as  four  percentage  points  in 
nozzle  efficiency  over  the  expected  range  of  maximum  Weber  number).  In 
contrast,  the  droplet  drag  coefficient  and  the  droplet  heat  transfer  co¬ 
efficient  were  found  to  have  a  lesser  influence  on  the  nozzle  performance. 
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7.  Performance  criteria  of  two-pha6e  nozzles  applicable  to  work  extraction 
systems  were  identified  to  be  nozzle  efficiency  and  mixture  velocity  at 
the  nozzle  exit.  But  the  performance  and  design  of  the  selected  nozzle 
must  be  carefully  integrated  with  the  overall  work  extraction  system  to 
yield  the  optimum  overall  system. 

8.  Using  a  simplified  sizing  procedure  based  on  isothermal  theory,  estimates 
of  off-design  gas  and  liquid  flow  rates  indicate  that  the  kinetic  energy 
profiles  at  off-design  can  be  made  relatively  flat  by  tailoring  the  flow 
rates  . 


R81-955229-4 


Two-Phase  Nozzle  Theory  and  Parametric  Analysis 
Phase  II  -  Parametric  Analysis  and  Optimisation 


The  overall  objective  of  this  analytical  study  program  is  to  advance  the 
basic  understanding  of  two-phase  nozzles,  particularly  from  the  energy  conversion 
point  of  view.  To  meet  this  overall  objective,  the  specific  objectives  of  the  Phase 
11  work  described  herein  are  to  conduct  parametric  analyses  in  terms  of  the  effect  of 
the  governing  parameters  on  two-phase  nozzle  performance  characteristics  and  to 
identify  optimum  design  criteria  for  two-phase  nozzles.  As  part  of  the  earlier  Phase 
1  work  (which  is  described  in  detail  in  Ref.  1),  a  working  theory  of  two-component 
two-phase  (gas-liquid)  nozzles,  identified  as  the  Model  with  Droplet  Heat  Transfer, 
was  developed  which  includes  droplet  breakup  and  the  interphase  transport  of  momentum 
and  thermal  energy  but  excludes  mass  transfer  and  frictional  losses.  This  working 
theory  was  computerized  in  a  form  suitable  for  parameteric  and  optimization  analyses 
of  these  two-phase  nozzles.  The  Phase  II  work  described  herein  consists  of  results 
obtained  from  four  study  tasks:  II. 1)  Sensitivity  Analysis  of  Droplet  Characteristics 
II. 2)  Parametric  Performance  Analysis;  II. 3)  Optimum  Nozzle  Designs;  and  11.4) 
Assessment  of  Work  Extraction  Potential.  The  results  obtained  in  each  of  these  four 
tasks  are  presented  in  the  following  sections. 


II. 1  Sensitivity  Analysis  of  Droplet  Characteristics 

The  sensitivity  of  the  two-phase  nozzle  performance  predictions  (obtained  from 
the  working  theory  that  was  developed  during  Phase  I  of  this  program  and  is  presented 
in  Ref.  1)  to  assumptions  and  simplifying  correlations  concerning  droplet  characteris¬ 
tics  was  examined  by  varying  the  correlations  and  then  comparing  the  change  in  the 
predicted  performance.  To  make  numerical  estimates  of  nozzle  performance  characteris¬ 
tics,  correlations  of  experimental  data  are  required  in  three  areas:  droplet  breakup, 
the  droplet  drag  coefficient,  and  the  droplet  heat  transfer  coefficient.  Estimates  of 
the  sensitivities  of  the  nozzle  performance  predictions  to  these  three  droplet  charac¬ 
teristics  were  made  in  this  portion  of  the  study. 

The  sensitivity  of  the  predicted  nozzle  efficiency  to  the  assumptions  made  for 
droplet  breakup  was  found  to  be  significant.  In  contrast,  the  sensitivities  of  the 
predicted  nozzle  efficiency  to  the  assumptions  for  the  droplet  drag  coefficient  and 
the  droplet  heat  transfer  coefficient  were  found  to  be  low,  in  part  because  of  the 
relatively  high  heat  transfer  coefficient  (from  the  liquid  droplet  to  the  gas) 
inherent  with  the  very  small  droplets. 
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II. 1.1  Droplet  Breakup 


As  discussed  in  Ref.  1  and  elsewhere,  droplet  breakup  in  a  two-phase  mixture 
is  governed  primarily  by  the  ratio  of  two  forces:  the  aerodynamic  pressure  forces, 
and  the  surface  tension  forces.  The  ratio  of  these  two  forces  is  known  as 
the  Weber  Number,  which  is  defined  as: 


We 


p  V  D 
8  slip 

2o 


(II. 1) 


where  is  the  gas  density,  is  the  slip  velocity  (vs^p  =  -  V^), 

D  is  the  droplet  diameter,  and  o  is  the  surface  tension  of  the  liquid.  The 
subscript  "i"  refers  to  the  liquid  phase,  and  the  subscript  "g"  refers  to  the  gas 
phase . 


In  the  computerized  model,  the  critical  value  of  the  Weber  Number,  We)cr£t» 
would  limit  the  maximum  droplet  size,  that  can  exist  at  the  local  flow 

conditions  in  a  two-phase  flow: 


2o  We) 


Crit 


MAX 


P  V 

slip 
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The  generally  accepted  value  of  We)cr^  is  6.0  (Ref.  2),  within  a  factor  of  about 
2  (for  example,  see  Ref.  3).  However,  recent  experiments  by  Alger  (Ref.  4)  with 
steam-water  nozzles  suggest  that  the  appropriate  value  of  We)cr£t  should  be 
approximately  1.0.  Areas  of  high  acceleration  in  the  nozzle  could  produce  greater 
droplet  deformation  as  compared  to  nearly  spherical  droplets  in  less  accelerated 
regions. 

The  sensitivity  of  the  nozzle  performance  predictions  to  droplet  breakup  in 
terms  of  the  exit  velocities  and  temperatures  and  the  nozzle  efficiency  is  shown 
in  Table  II. 1  for  three  values  of  the  critical  Weber  Number.  These  values  generally 
cover  the  range  found  experimentally  by  previous  investigators  and  represent  the 
range  over  which  the  sensitivity  of  nozzle  performance  predictions  is  of  interest. 

In  Table  II. 1,  V  is  velocity,  r  is  the  loading  ratio  (r  =  mt/m^)  where 
m  is  mass  flow  rate,  T  is  temperature,  and  n  is  the  nozzle  efficiency  (which 
is  defined  later  in  Eq.  II. 9).  Droplet  breakup  is  seen  to  have  a  significant  effect 
on  the  predicted  efficiency  because  the  smaller  droplets  have  a  lower  slip  velocity 
(V  -Vj)  and  a  higher  rate  of  heat  transfer  to  the  gas  because  of  the  increased 
total  surface  area  of  the  droplets.  At  the  selected  operating  conditions  (given  in 
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Table  II. 1)  of  a  nozzle  with  a  two-phase  mixture  of  steam  and  Dovtherm  A  (DTA), 
decreasing  the  value  of  the  critical  Weber  Number  by  a  factor  of  six  results  in  an 
increase  in  the  predicted  nozzle  efficiency  of  2.7  percentage  points  at  a  loading 
ratio  of  one,  while  at  a  loading  ratio  of  AO,  the  increase  in  predicted  nozzle 
efficiency  is  4.3  percentage  points. 

11.1.2  Droplet  Drag 

As  discussed  in  Ref.  I,  the  empirical  data  base  for  the  drag  coefficient 
(Cp)  was  developed  from  previous  investigations  published  in  the  literature,  and 
the  comprehensive  expression  that  was  selected  for  use  in  the  computerized  model 
is : 


*n  c:  *=  3.271  -  0.8893  L  +  0.03417  L*  +  0.001443 
D  Re  Re  R 


where  LRe  =  In  (ReRej).  ReReJ  is  the  relative  Reynolds  number 

(ReRel  2  pg  Vslip  D/,pR)’  and  pg  is  the  838  viscosity*  Stonecypher 

(as  cited  in  Ref.  5)  developed  this  expression  by  making  a  least-squares  fit  of 

the  drag  data  for  spherical  particles  that  was  tabulated  by  Perry,  et  al.  (Ref.  6). 

Using  the  computerized  model,  the  sensitivity  of  the  nozzle  performance 
predictions  to  the  droplet  drag  coefficient  was  calculated  and  is  shown  in  Table 
II. 2  for  three  values  of  the  coefficient:  1)  80  percent  of  the  value  predicted  by 
the  drag  coefficient  correlation  of  Eq .  II. 2  at  the  local  relative  Reynolds  number 
for  each  of  the  axial  nozzle  segments,  2)  100  percent  of  this  value,  and  3)  120 
percent  of  this  value.  This  range  of  values  (plus  and  minus  20  percent  of  the 
nominal  value)  brackets  a  reasonable  range  of  scatter  that  might  exist  in  the 
original  experimental  drag  data.  The  sensitivity  of  the  calculated  exit  velocities 
and  temperatures  was  found  to  be  small  at  a  loading  ratio  of  40,  with  only  a  0.7 
percentage  point  increase  in  the  calculated  efficiency  as  the  drag  coefficient  is 
increased  over  the  range.  Conversely,  at  a  loading  ratio  of  one,  the  sensitivity 
of  the  calculated  nozzle  efficiency  was  found  to  be  extremely  low  because  the 
relative  impact  of  the  liquid  is  less,  and  the  sensitivities  of  the  exit  velocities 
and  temperatures  are  greater  than  at  higher  values  of  the  loading  ratio  but  are 
still  relatively  low. 

1 1 . 1 . 3  Droplet  Heat  Transfer 

The  heat  transfer  coefficient  (h)  between  the  liquid  droplets  and  the  gas 
is  numerically  evaluated  in  the  computerized  model,  as  discussed  in  Ref.  1,  with  a 
dimensionless  correlation  recommended  by  McAdams  (Ref.  7)  for  relative  Reynolds 
numbers  in  the  range  of  25  to  10,000: 
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where  is  the  thermal  conductivity  of  the  gas.  For  relative  Reynolds  numbers 
in  the  range  of  1.0  to  25,  the  expression  recommended  by  Kreith  (Ref.  8)  is  used 
to  calculate  the  heat  transfer  coefficient: 
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where  C  „  is 
Pg 


the  specific  heat 


(at  constant  pressure)  of  the  gas. 


Using  the  computerized  model,  the  sensitivity  of  the  nozzle  performance 
predictions  to  the  heat  transfer  coefficient  was  calculated  and  is  shown  in  Table 
II. 3  for  three  values  of  the  coefficient:  1)  80  percent  of  the  value  predicted  by 
the  heat  transfer  coefficient  correlations  of  Eqs.  II. 3  and  II. 4  at  the  local 
relative  Reynolds  number  for  each  of  the  axial  nozzle  segments,  2)  100  percent  of 
this  value,  and  3)  120  percent  of  this  value.  This  range  of  values  brackets  the 
♦  20-percent  scatter  that  is  typical  of  heat  transfer  data.  The  sensitivity 
of  Che  calculated  efficiency  and  the  exit  velocities  and  temperatures  to  heat 
transfer  was  found  to  be  low  over  the  loading  ratio  range  of  1.0  to  40.,  and 
decreases  as  the  loading  ratio  is  increased.  The  low  sensitivity  is  attributable, 
in  part,  to  the  fact  that  the  heat  transfer  coefficient  is  already  very  high 
because  the  liquid  droplets  are  very  small  in  diameter  at  the  operating  conditions 
of  Table  II. 3. 
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II. 2  Parametric  Performance  Analysis 

Using  the  working  theory  (identified  as  the  Model  with  Droplet  Heat  Transfer) 
that  was  previously  developed  on  Phase  I  of  this  program  (Ref.  1),  parametric 
performance  analyses  were  made  for  selected  two-phase  nozzles  in  terms  of  the 
operating  conditions.  The  effects  of  all  the  governing  parameters  on  nozzle 
performance  were  investigated,  and  efficiency  maps  of  a  two-phase  nozzle  were 
constructed  based  on  these  analyses.  The  power  split  between  the  liquid  and  the 
gas  streams  was  estimated  under  different  operating  conditions.  Table  II. 4 
summarizes  the  relative  effect,  qualitatively  rated,  of  the  governing  flow  param¬ 
eters  on  nozzle  performance  predictions  in  terms  of  the  exit  temperatures  and 
velocities  of  each  phase  and  of  the  nozzle  efficiency.  The  parameters  that  have 
the  most  significant  effect  on  nozzle  efficiency  are:  loading  ratio,  nozzle 
pressure  ratio  (or  exit  pressure),  gas  density,  liquid  surface  tension  and  the 
related  droplet  size,  and  nozzle  length.  The  results  of  the  parametric  analyses 
which  form  the  base  for  the  ratings  of  the  relative  effects  shown  in  Table  11.4 
are  presented  and  discussed  in  this  section.  Generally,  the  value  of  the  Weber 
Number  was  set  at  6.0,  the  value  of  the  inlet  gas  velocity  was  set  at  100  ft/sec, 
the  value  of  the  inlet  liquid  velocity  was  set  at  70  ft/sec,  and  a  nozzle  length 
of  six  inches  was  selected  for  the  parametric  analyses,  unless  otherwise  stated. 

II. 2.1  Effect  of  Liquid/Gas  Ratio 

The  effect  of  the  liquid/gas  ratio  was  investigated  in  terms  of  the  ratios  of 
mass,  volume,  and  thermal  capacity  of  the  two  streams.  The  mass  ratio  of  the  two 
streams  is  generally  known  as  the  loading  ratio  (r  s  rn^/m  ) ,  while  the 
ratios  of  volume  and  thermal  capacity  are  directly  related  to  this  loading  ratio. 
For  two-component  two-phase  nozzle  flow  without  mass  transfer  between  the  two 
phases  (which  is  the  situation  for  the  working  theory  developed  during  Phase  I  of 
this  program),  the  loading  ratio  is  constant  throughout  the  nozzle,  but  the 
ratios  of  volume  and  thermal  capacity  will  vary  as  the  two-phase  mixture  flows 
through  the  nozzle,  depending  on  the  effect  of  local  temperature  and  pressure  on 
the  density  and  specific  heat  of  the  gas  and  the  liquid.  The  volume  ratio  is 
defined  as  the  ratio  of  the  volumetric  flows  of  the  two  streams: 


r 

v 


m  /c 
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(II. 5) 


where  ot  is  the  density  of  the  liquid.  The  thermal  capacity  ratio  is  corres¬ 
pondingly  defined  as  the  ratio  of  the  thermal  capacities  (per  unit  temperature 
rise)  of  the  two  streams: 
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where  C  is  the  specific  heat  of  the  liquid.  This  definition  is  analogous  to  the 
definition  of  capacity  ratio  in  heat  exchanger  design  theory  (Ref.  9). 

Figure  II. 2  shows  the  effect  of  the  (mass)  loading  ratio  on  performance 
predictions  of  nozzle  throat  and  exit  conditions.  P  is  the  local  pressure,  and 
the  subscript  "in"  refers  to  the  nozzle  inlet  while  the  subscript  "exit"  refers  to 
the  nozzle  exit.  As  the  loading  ratio  is  reduced,  both  throat  and  exit  velocities 
of  both  phases  increase,  and  both  throat  and  exit  temperatures  of  both  phases 
decrease.  These  situations  are  explained  by  reasoning  that  the  effect  of  the  liquid 
phase  decreases  as  the  loading  ratio  is  decreased,  and  hence  the  gas  phase  tends 
to  behave  increasingly  as  a  pure  gas.  The  slip  velocity  increases  as  the  loading 
ratio  is  reduced  and  the  droplet  size  decreases  as  a  result,  because  the  higher 
6lip  velocity  means  higher  aerodynamic  forces  which  tend  to  break  up  the  droplets. 
The  nozzle  efficiency  goes  through  a  minimum,  at  first  decreasing  because  of  the 
higher  slip  when  the  loading  ratio  is  reduced,  but  then  increasing  because  the 
presence  of  the  liquid  phase  has  a  lesser  effect  when  the  loading  ratio  is  reduced 
further . 

As  indicated  above,  both  the  volume  ratio  and  the  thermal  capacity  ratio  are 
proportional  to  the  mass  loading  ratio.  Hence,  the  effects  of  the  volume  and 
thermal  capacity  ratios,  evaluated  at  the  nozzle  inlet  conditions  (of  temperature 
and  pressure)  are  the  same  as  the  effect  of  the  mass  loading  ratio,  and  are  shown 
in  Fig.  II. 2. 

The  exit  specific  power  is  shown  in  Fig.  II. 3  as  a  function  of  (mass)  loading 
ratio  and  pressure  ratio  (PR  ■  pin^exit^’  The  ex^t  sPecific  power  is  the 
time  rate  of  flow  of  kinetic  energy  of  the  two-phase  mixture  at  the  nozzle  exit, 
divided  by  the  exit  flow  area  (Afix^t): 


Specific  Power) 


1  .  -2 
-  m  V 
2  exit 


(II. 7) 


where  m  is  the  total  mass  flow  rate  (m  ■  m^/  +  m^)  and  ^exit 

is  the  mass-weighted  exit  velocity.  The  mass-weighted  velocity  is  defined  as: 

V  *  (V  ♦  r  Vt)/(1  +  r).  As  seen  in  Fig.  II. 3,  the  specific  power  increases 
when  either  the  loading  ratio  or  the  pressure  ratio  is  decreased.  As  the  loading 
ratio  is  decreased,  the  required  exit  area  decreases  because  the  gas  exit  velocity 
is  increasing  (as  seen  in  Fig.  II. 2),  and  less  flow  area  is  needed  to  pass  the 
higher-velocity  mixture;  with  less  liquid  relative  to  the  gas  mass,  the  gas  exit 
velocity  is  higher  at  the  lower  loading  ratio.  As  the  pressure  ratio  decreases 
(and  hence  the  exit  pressure  increases,  for  fixed  inlet  pressure),  the  exit  gas  is 
denser  and  less  exit  area  is  needed  to  pass  the  same  mass  flow  at  a  given  loading 
ratio.  With  less  nozzle  exit  flow  area,  then,  the  specific  power  increases. 
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II. 2. 2  Pressure  and  Temperature  Effects 

Pressure  and  temperature  effects  on  nozzle  performance  predictions  were 
investigated  parametrically  in  terms  of  inlet  pressure,  exit  pressure,  pressure 
ratio  (P^n/Pex^t) ,  inlet  temperature,  and  inlet  temperature  difference  between 
the  gas  and  the  liquid. 

Pressure  Effects 


With  the  exit  pressure  held  constant.  Fig.  II. 4  shows  the  effect  of  inlet 
pressure  on  nozzle  performance  predictions  for  a  two-phase  mixture  of  steam  and 
DTA,  for  two  values  of  loading  ratio  which  span  the  range  of  interest  for  propulsion 
applications.  The  effect  on  nozzle  efficiency  is  small,  with  the  nozzle  efficiency 
increasing  slightly  as  the  inlet  pressure  is  increased.  The  exit  velocities 
increase  and  the  exit  temperatures  decrease  as  the  inlet  pressure  is  increased, 
because  the  gas  phase  (steam)  is  being  expanded  through  a  larger  pressure  ratio. 

Figure  II. 5  shows  the  parametric  effect  of  nozzle  exit  pressure,  with  a  fixed 
inlet  pressure  of  1000  psia.  As  the  exit  pressure  is  decreased,  the  gas  phase  is 
expanded  more,  and  the  gas  exit  velocity  increases  substantially  while  the  liquid 
exit  velocity  increases  by  a  lesser  amount.  Hence,  the  slip  velocity  (V^-Vj) 
between  the  two  phases  increases,  and  the  nozzle  efficiency  decreases  as  a  result. 

In  addition,  with  the  greater  expansion  of  the  gas  phase  when  the  exit  pressure  is 
decreased,  the  gas  exit  temperature  also  decreases  substantially.  Here,  the 
liquid  phase  acts  as  a  heat  source  to  the  gas  phase,  and  heat  transfer  from  the 
hotter  liquid  retards  this  gas  temperature  decrease,  but  the  short  residence  time 
of  the  two-phase  mixture  in  the  nozzle  does  limit  the  amount  of  heat  that  can  be 
transferred  from  the  liquid  to  the  gas.  At  the  lower  loading  ratio  (r  *  10),  the 
liquid  exit  temperature  is  lower  and  the  effect  of  the  lower  exit  pressure  is 
larger,  in  both  instances  because  the  liquid  heat  source  is  relatively  smaller 
than  at  a  loading  ratio  of  40. 

Shown  in  Fig.  II. 6  is  the  effect  of  pressure  ratio  (Pjn/Pexit)»  over  an 
extremely  wide  range,  at  a  fixed  inlet  pressure.  As  the  pressure  ratio  is  increased, 
both  the  efficiency  and  the  exit  gas  temperature  decrease,  the  exit  gas  velocity 
rises  sharply,  while  the  exit  liquid  temperature  falls  and  the  exit  liquid  velocity 
rises  only  for  pressure  ratios  less  than  about  50  and  are  otherwise  little  affected 
by  the  pressure  ratio.  For  the  results  of  Fig.  11.6,  the  pressure  ratio  is 
increased  by  decreasing  the  exit  pressure  to  lower  values.  For  example,  at  a 
pressure  ratio  of  400  and  an  inlet  pressure  of  600  psia,  the  exit  pressure  is  1.5 
psia.  Figure  II. 7  shows  a  nozzle  design  at  this  low  exit  pressure.  In  the  last 
0.2  ft  of  the  nozzle,  the  gas  temperature  drops  rapidly  while  the  liquid  temperature 
drops  only  moderately,  and  the  gas  velocity  rises  sharply  while  the  liquid  velocity 
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rises  only  moderately.  In  general,  at  low  exit  pressures,  for  a  pressure  profile 
linear  with  nozzle  length,  the  local  expansion  ratio  near  the  exit  is  large 
relative  to  design  conditions  that  have  higher  exit  pressures.  This  large  local 
expansion  ratio  then  results  in  large  increases  in  gas  velocity  and  large  decreases 
in  gas  temperature  as  the  gas  expands.  The  residence  time  is  too  short  at  the 
high  velocities  (both  for  momentum  transfer  from  the  gas  to  the  liquid  and  for 
heat  transfer  from  the  liquid  to  the  gas),  to  narrow  the  velocity  and  temperature 
differences  between  the  gas  and  the  liquid,  respectively.  The  effect  of  exit 
pressure  on  nozzle  shape  is  shown  in  Fig.  II. 8,  in  terms  of  the  relative  nozzle 
diameter  at  the  throat  and  at  the  exit,  for  a  six-inch-long  nozzle  with  the 
pressure  profile  linear  with  length.  Exit  pressure  has  no  effect  on  throat 
diameter,  but  the  exit  diameter  must  be  increased  to  accommodate  the  lower-density 
exit  gas  at  the  lower  exit  pressure,  even  though  the  exit  gas  velocity  increases 
as  the  exit  pressure  is  reduced;  this  trend  of  gas  velocity  was  seen  in  Fig. 

II. 6.  In  general,  the  nozzle  exit  area  (which  is  proportional  to  dgX^t)  in 
Fig.  II. 8  is  increasing  faster  than  the  inverse  of  the  exit  gas  density  (this 
inverse  is  proportional  to  1/Pex£t)  because  the  exit  gas  velocity  is  not  increas¬ 
ing  as  fast  as  the  exit  gas  density  is  decreasing. 

Temper ature  EfXefLc_5. 

Table  II. 5  shows  the  effect  of  inlet  temperature  on  nozzle  performance 
predictions,  at  two  values  of  the  loading  ratio.  At  a  given  loading  ratio,  the 
effect  of  inlet  temperature  on  the  liquid  temperature  drop  through  the  nozzle 
[Tt)£n  -  Tt)ex^t]  is  6mall,  as  is  the  effect  on  the  temperature  differential 
between  the  liquid  and  the  gas  at  the  nozzle  exit.  As  the  inlet  temperature  is 
reduced  from  1160  R  to  960  R,  the  liquid  and  gas  velocities  and  the  nozzle  efficiency 
decrease,  while  the  slip  velocity  increases  because  the  liquid  droplets  are  larger 
at  the  lower  temperature  where  the  liquid  surface  tension  is  larger.  With  the 
higher  slip  velocity,  the  nozzle  efficiency  is  expected  to  be  lower. 

The  effect  of  the  inlet  temperature  difference  on  nozzle  performance  predic¬ 
tions  is  shown  in  Fig.  II. 9.  The  only  major  effect  of  this  temperature  difference 
is  seen  in  the  exit  temperatures  of  the  gas  and  the  liquid,  where  the  exit  tempera¬ 
tures  fall  when  the  gas  temperature  is  reduced  to  provide  the  temperature  difference. 

This  trend  is  expected  from  enthalpy  considerations.  The  nozzle  efficiency  also 
decreases  slightly  as  the  inlet  temperature  is  increased,  because  the  colder 
two-phase  mixture  means  that  the  liquid  surface  tension  is  higher  and  hence  the 
liquid  droplets  are  larger  according  to  the  Weber  Number  breakup  criterion  (Eq. 

II. la).  This  data  on  the  effect  of  inlet  temperature  difference  is  replotted  in 
Fig.  II. 10  as  a  function  of  loading  ratio.  The  effect  is  relatively  larger  at  the 
lower  values  of  the  loading  ratio,  where  the  gas  phase  is  a  larger  fraction  of  the 
two-phase  mixture. 
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II. 2, 3  Effects  of  Theraophysical  Properties 

The  effects  of  theraophysical  properties  of  the  gas  and  of  the  liquid  on 
nozzle  performance  characteristics  were  investigated  by  judiciously  varying  these 
properties . 

Ef_fec t s^  £f_Gj»s_Theraophy£ica_l  Properties 

Table  II. 6  shows  the  effect  of  the  gas  specific  heat,  at  three  values  relative 
to  the  specific  heat  of  steam:  twice  the  specific  heat  of  steam,  equal  to  that, 
and  half  of  that  value.  As  expected,  the  effect  on  exit  velocities,  which  are 
primarily  governed  by  momentum  transfer,  is  negligible.  The  effect  on  nozzle 
efficiency  is  small  and  decreasing  as  the  loading  ratio  is  increased,  and  the 
effect  on  exit  temperatures  is  also  small,  in  part  because  the  gas  phase  is  only 
one-eleventh  of  the  total  mass  at  a  loading  ratio  of  ten,  for  example. 

The  effect  of  gas-phase  density  on  nozzle  performance  predictions  is  shown  in 
Fig.  11.11.  The  calculations  were  made  by  varying  the  molecular  weight  of  the  gas 
(hence  the  gas  density  which  is  calculated  with  the  perfect  gas  law),  while  the 
relationships  for  the  other  theraophysical  properties  (viscosity,  thermal  conduc¬ 
tivity,  and  specific  heat)  of  the  gas  were  not  changed  and  were  those  of  steam. 

As  the  gas  density  increases,  the  nozzle  efficiency  increases,  and  the  exit 
velocities  and  the  slip  velocity  (Vg-Vt)  both  decrease.  The  lower  velocities  mean 
both  a  longer  residence  time  in  the  nozzle  for  increased  heat  transfer  from 
the  liquid  to  the  gas,  and  higher  exit  temperatures  for  both  phases  because  a 
lower  fraction  of  the  constant  total  energy  is  in  the  form  of  kinetic  energy.  If 
the  other  theraophysical  properties  of  the  gas  remain  unchanged,  a  heavier  gas, 
therefore,  results  in  higher  nozzle  efficiencies. 

Ef_fect£  of_L_i_quid[  Thermo s_i c £  1_P T^ope rt _ie 

The  effects  of  the  theraophysical  properties  of  the  liquid  were  also  investi¬ 
gated,  and  Fig.  11.12  shows  the  relative  effect  of  liquid  specific  heat  on  nozzle 
performance  characteristics  (all  other  liquid  properties  being  those  of  Dowtherm 
A).  Three  values  of  liquid  specific  heat  were  investigated:  the  specific  heat 
of  Dowtherm  A,  and  half  and  twice  that  value.  At  a  loading  ratio  of  ten,  halving 
the  liquid  specific  heat  results  in  only  a  0.15  percentage  increase  in  the 
predicted  nozzle  efficiency,  and  the  effect  is  smaller  at  larger  values  of  the 
loading  ratio.  Of  course,  the  effect  on  the  calculated  exit  temperatures  of  the 
liquid  and  of  the  gas  are  much  larger:  the  liquid  with  the  lower  specific  heat 
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has  a  lower  exit  temperature  and  the  accompanying  gas  exit  temperature  is  also 
lower,  because  the  liquid  is  a  relatively  smaller  heat  source  at  the  same  loading 
ratio  than  is  a  liquid  with  a  higher  specific  heat.  Overall,  the  effect  of  liquid 
specific  heat  on  the  nozzle  efficiency,  hence  its  performance,  is  small. 

Figure  11.13  shows  the  effect  of  liquid  Burface  tension  on  nozzle  performance 
predictions.  As  discussed  previously  in  connection  with  the  Weber  Number  breakup 
criterion  (Eq.  II. la),  liquid  surface  tension  is  a  primary  factor  in  the  determina¬ 
tion  of  the  droplet  size.  And,  as  the  surface  tension  is  decreased,  the  liquid 
droplets  will  break  up  into  smaller  sizes.  With  smaller  droplets,  at  a  given 
loading  ratio,  the  velocity  slip  (Vg-V^)  between  the  two  phases  is  lower,  the  heat 
transfer  between  the  liquid  and  the  gas  is  higher  because  there  i6  more  droplet 
surface  area,  and  the  nozzle  efficiency  is  higher  as  a  result.  These  are  the 
trends  displayed  as  a  function  of  loading  ratio,  in  Fig.  11.13,  which  shows 
results  for  three  values  of  the  liquid  surface  tension:  the  real  value  of  DTA, 
half  of  that  value,  and  twice  that  value.  Except  for  the  surface  tension,  the 
thermophysical  properties  of  the  liquid  phase  are  those  of  DTA,  and  the  thermophysical 
properties  of  the  gas  phase  are  those  of  steam.  If  the  other  thermophysical 
properties  of  the  liquid  remain  unchanged,  lower  liquid  surface  tension,  therefore, 
results  in  higher  nozzle  efficiency. 

The  effect  of  liquid  density  is  shown  in  Fig.  11.14  in  terms  of  three  values 
of  density:  the  density  of  DTA,  and  both  20  percent  above  and  20  percent  below 
the  value  for  the  DTA  density.  All  other  liquid  properties  are  those  of  DTA.  The 
effect  of  liquid  density  on  the  gas  and  liquid  exit  temperatures  is  negligible  and 
is  not  shown.  At  a  given  value  of  the  loading  ratio,  the  lower  liquid  density 
leads  to  higher  velocities,  lower  slip  velocity  (V  -V^),  and  higher  nozzle  effi¬ 
ciency.  With  a  given  amount  of  momentum  transfer  from  the  gas  to  the  liquid, 
the  lower  liquid  density  means  that  the  liquid  velocity  will  be  higher.  For  a  20 

percent  decrease  in  liquid  density,  the  nozzle  efficiency  increases  by  1.6  percent¬ 

age  points  at  a  loading  ratio  of  40,  and  by  1 . 1  percentage  points  at  a  loading 
ratio  of  10. 

II. 2. 4  Effects  of  Inlet  Velocity,  Droplet  Size,  and  Nozzle  Length 

The  effect  of  inlet  velocity  on  nozzle  performance  predictions  is  summarized 
in  Table  II. 7  and  is  seen  to  be  small.  Part  a)  of  this  table  presents  nozzle  exit 
results  at  three  levels  of  inlet  gas  velocity,  with  the  velocity  ratio  (Vt/Vg)  held 

constant  at  0.7,  for  three  values  of  the  loading  ratio.  Part  b)  shows  results 

for  a  loading  ratio  of  20,  but  at  three  values  of  the  inlet  velocity  ratio,  0.7, 

0.8,  and  0.9,  at  two  levels  of  inlet  gas  velocity.  In  all  cases,  the  effect  on 
exit  temperatures  is  negligible.  As  the  inlet  velocity  is  increased,  the  nozzle 
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efficiency  increases  slightly,  as  do  the  exit  velocities,  but  by  much  smaller 
amounts  than  the  increase  in  inlet  velocity.  At  a  given  inlet  gas  velocity,  the 
effect  of  velocity  ratio  is  negligible.  Hence,  in  general,  the  effect  of  inlet 
velocity  on  nozzle  performance  predictions  is  minor. 

The  effect  of  droplet  size  is  shown  in  Fig.  11.15.  As  expected  from  discus¬ 
sions  above,  the  effect  of  droplet  size  on  nozzle  efficiency  is  significant 
because  of  its  impact  on  the  transfer  of  momentum  and  heat.  The  results  in  Fig. 

11.15  are  presented  in  terms  of  the  maximum  Weber  Number,  because  this 

is  the  factor  which  is  postulated  to  control  droplet  size  in  the  Model  with 
Droplet  Heat  Transfer.  The  droplets  break  up  until  their  diameter  results  in  a 
droplet  Weber  Number  equal  to  the  maximum.  A  smaller  value  of  the  maximum  Weber 
Number  results  in  smaller  droplet  sizes,  as  indicated  by  the  results  in  Fig.  11.15 
for  droplet  diameter  (D),  but  the  droplet  diameter  does  not  vary  proportionately 
with  the  maximum  Weber  Number  because  the  slip  velocity  also  changes  when  the 
droplet  size  changes.  The  droplet  drag  force,  which  is  influenced  by  the  droplet 
size,  couples  the  gas  to  the  liquid  and  is  indirectly  related,  through  the  force 
balance  on  the  droplet,  to  the  velocity  difference  between  the  two  phases.  As 
seen  in  Fig.  11.15,  the  major  effect  of  the  droplet  size  is  to  increase  the  nozzle 
efficiency  by  about  two  percentage  points  when  the  critical  Weber  Number  (for 
droplet  breakup)  is  cut  by  a  factor  of  two,  and  the  corresponding  droplet  size  is 
then  about  25  percent  smaller.  The  smaller  droplets  result  in  a  higher  velocity 
ratio  (V  / V, ) ,  which  means  lower  slip  velocity  and  higher  nozzle  efficiency. 

Figure  11.16  shows  the  effect  of  nozzle  length  on  performance  predictions. 

The  main  effect  is  on  nozzle  efficiency  which  increases  as  nozzle  length  increases 
because  the  residence  time  of  the  two-phase  mixture  is  longer.  And  the  longer 
residence  time  means  more  transfer  of  momentum  and  heat  between  the  two  phases. 

Hence,  as  the  nozzle  length  is  increased,  the  velocity  slip  (Vg-V^)  is  reduced,  as 
is  the  temperature  difference  between  the  hot  liquid  and  the  cold  gas.  As  mentioned 
previously,  the  Model  with  Droplet  Heat  Transfer  neglects  friction,  and  this  is  a 
reasonable  assumption  for  shorter  nozzles  to  be  used  in  propulsion  applications  where 
space  is  likely  to  be  limited. 

11.2.5  Power  Split  of  Two-Phase  Nozzles 


The  power  split  of  the  nozzle  is  defined  as  the  ratio  of  the  kinetic  energy 
of  the  liquid  fraction  at  the  nozzle  exit  to  the  kinetic  energy  of  the  gas  fraction 
at  the  nozzle  exit: 
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where  K  • r  is  the  exit  value  of  the  velocity  ratio  (K  =  V  /V  ) .  Because 
of  the  assumption  that  no  mass  transfer  occurs  between  the  gas  and  the  liquid 
phases,  the  mass  loading  ratio  does  not  change  from  the  inlet  value.  Therefore, 
to  estimate  the  power  split,  the  inlet  loading  ratio  and  the  velocity  ratio  at  the 
nozzle  exit  must  be  obtained.  The  drag  force  on  the  liquid  droplets  accelerate 
the  droplets,  but  the  droplet  velocity  lags  the  gas  velocity  throughout  the 
nozzle.  The  velocity  ratio  at  the  nozzle  exit  is  calculated  as  one  of  the  results 
from  the  computerized  Model  with  Droplet  Heat  Transfer  and  was  shown  earlier  in 
the  parametric  results  (as  in  Fig.  11.2).  Generally,  the  range  of  values  of 
Kexit  to  about  0.5  to  0.9,  but  the  value  of  the  loading  ratio  i6  much 

higher,  ranging  from  zero  (all-gas  flow)  to  50  and  higher,  depending  on  the  design 
conditions.  Hence,  the  loading  ratio  which  is  one  part  of  the  specified  inlet 
operating  conditions  is  the  more  important  factor  in  the  estimation  of  the  power 
split.  At  a  loading  ratio  of  40  for  the  steam/DTA  mixture  and  the  nozzle  inlet 
conditions  specified  in  Fig.  11.2,  the  value  of  the  exit  velocity  ratio  i6  Kex^t  *  0.78 
and  the  power  split  is  then  24.4,  as  estimated  with  Eq.  11.8.  As  expected,  the 
bulk  of  the  kinetic  energy  is  carried  at  this  loading  ratio  by  the  liquid  phase 
despite  the  higher  exit  velocity  of  the  gas  phase.  If  the  pressure  ratio  is  increased 
to  40  with  the  loading  ratio  still  at  40,  then  the  exit  velocity  is  about  0.72  and  the 
power  split  is  20.7.  Figure  11.17  shows  the  power  split  as  a  function  of  the  loading 
ratio.  Also  shown  in  Fig.  11.17  is  the  same  information  replotted  as  the  liquid 
energy  fraction,  which  is  the  fraction  of  the  total  kinetic  energy  of  the  mixture  that 
is  carried  by  the  liquid  phase. 


II. 2. 6  Efficiency  Maps  of  Two-Phase  Nozzles 


Based  on  the  parametric  performance  analyses  discussed  above,  efficiency  maps 
for  a  nozzle  with  a  two-phase  mixture  of  steam  and  Dowtherm  A  (DTA)  were  constructed. 
As  noted  previously  (Ref.  1),  the  isentropic  efficiency  is  a  measure  of  nozzle 
efficiency  and  is  defined  as: 


1  i  (V~  -V  )  V  -V, 

2  exit  inlet  exit  Inlet  (II. 9) 

71  m  Ah)  2  Ah) 

max  max 

This  efficiency  is  the  ratio  of  the  change  in  kinetic  energy  per  unit  mass  of  the 
mixture  to  the  ideal  change  of  enthalpy  (into  kinetic  energy)  for  the  isentropic 
expansion  of  a  homogeneous  two-phase  mixture.  (For  operating  conditions  where  the 
nozzle  inlet  velocity  )  is  very  low  when  compared  to  the  exit  velocity  (Vex£t), 

then  ^exjt  >>  ^inlet  an<*  t*ie  ^n^et  velocity  term  can  be  neglected.) 


The  maximum  enthalpy  change  through  the  nozzle,  A^)max»  is  given  by  the 
expression: 
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where  AP  is  the  pressure  drop  through  the  nozzle.  Efficiency  maps  are  shown  in 
Figs.  II. 18,  11.19,  and  11.20  in  terms  of  pressure  effects,  temperature  effects, 
and  inlet  temperature  difference  effects,  respectively. 

Ef^fec t_of_  j^resjsure_Rat_io 

The  effect  of  pressure  ratio  (PR  *  Pin/Pexit)  on  nozzle  efficiency  at  two 
levels  of  inlet  pressure  is  seen  in  Fig.  11.18  over  a  range  of  loading  ratio.  At 
a  given  inlet  pressure,  as  the  exit  pressure  is  reduced  (and  hence  PR  increased), 
the  nozzle  efficiency  decreases.  At  the  lower  exit  pressure,  the  gas  expands  to  a 
higher  exit  velocity,  but  the  liquid  velocity  does  not  increase  by  a  proportionate 
amount  because  momentum  transfer  from  the  gas  to  the  liquid  is  a  rate  process 
which  therefore  requires  time  to  reach  equilibrium.  Hence,  at  the  lower  exit 
pressure,  because  the  residence  time  of  the  mixture  in  the  nozzle  is  shorter,  less 
momentum  is  transferred  from  the  gas  to  the  liquid  and  the  nozzle  efficiency  is 
lower.  At  a  given  pressure  ratio  but  at  a  lower  inlet  pressure,  the  exit  pressure 
is  lower,  and  the  same  rationale  of  less  momentum  transfer  also  explains  the  trend 
to  lower  efficiency. 

Effect  of  Inlet  Temperature 


In  Fig.  11.19,  the  increase  in  efficiency  as  the  inlet  temperature  is  increased 
can  be  explained  by  the  smaller  droplet  sizes  at  the  higher  temperature.  The 
liquid  surface  tension  (which  in  part  determines  the  maximum  droplet  size  through 
the  critical  Weber  Number  for  droplet  breakup,  Eq.  II. la)  decreases  as  the  tempera¬ 
ture  increases.  Hence,  the  droplets  are  smaller,  momentum  transfer  between  the 
smaller  droplets  and  the  gas  is  better,  the  exit  liquid  velocity  is  therefore 
higher,  and  the  nozzle  efficiency  is  higher  as  a  result. 

Effect _of^  InJ_e Temperature  ,Di_f  fejenc_e 

Figure  11.20  shows  that  for  a  constant  liquid  inlet  temperature,  the  effect 
of  inlet  temperature  difference  on  nozzle  efficiency  is  small  for  values  of  the 
loading  ratio  greater  than  about  five,  above  which  point  the  liquid  phase  dominates 
the  exit  temperature  of  the  two-phase  mixture.  Near  the  nozzle  inlet,  the  inlet 
temperature  difference  will  have  a  greater  effect,  however,  on  local  conditions, 
as  seen  earlier  in  the  Phase  I  work  (Fig.  1.26  of  Ref.  1).  Thus,  the  inlet 
temperature  itself  has  a  much  larger  effect  on  nozzle  efficiency  than  the  inlet 
temperature  difference. 

II. 2. 7  Comparisons  with  Experimental  Results 

As  discussed  in  Ref.  1,  experimental  data  have  been  taken  for  a  number  of 
two-phase  nozzles,  but  only  some  of  this  data  was  taken  with  a  two-component 
mixture  -  the  type  of  mixture  to  which  thei  Model  with  Droplet  Heat  Transfer 
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applies.  Figures  11.21,  11.22  and  11.23  compare  calculated  results  from  the  Model 
with  Droplet  Heat  Transfer  to  experimental  measurements  made  by  Elliott  and 
Weinberg  (Ref.  5)  at  three  levels  of  inlet  pressure  with  a  water-nitrogen  mixture 
flowing  through  a  50-inch-long  nozzle.  The  Model  with  Droplet  Heat  Transfer  is  a 
design  model,  and  as  such  is  normally  used  to  plan  families  of  nozzles  for  one  set 
of  operating  conditions.  Hence,  the  results  from  the  Model  with  Droplet  Heat 
Transfer  in  these  three  figures  are  actually  a  series  of  nozzle  designs,  over  a 
range  of  loading  ratios;  the  profile  of  pressure  drop  was  selected  to  be  linear 
with  axial  distance  through  the  nozzle.  Also  shown  in  these  three  figures  are 
performance  lines  predicted  by  Elliott  and  Weinberg  with  their  models  (Ref.  5)  for 
the  actual  nozzle  contour.  The  Elliott  and  Weinberg  prediction  labeled  "Real"  is 
based  on  a  model  that  includes  velocity  slip,  droplet  breakup,  and  heat  transfer 
but  no  friction;  while  their  prediction  labeled  "Friction"  includes  friction  in 
terms  of  wall  boundary  layer  losses. 

The  design-point  performances  shown  in  Figs.  11.21,  11.22  and  11.23  as 
calculated  for  a  series  of  nozzles  by  the  Model  with  Droplet  Heat  Transfer  compare 
well  with  the  experimental  data  of  Elliott  and  Weinberg  and  show  the  same  trend  as 
the  off-design  performance  of  one  nozzle. 

As  seen  by  the  definition  of  Eq.  11.9,  the  nozzle  efficiency  is  essentially  a 
function  of  the  mass-weighted  velocity  at  the  nozzle  exit.  Hence,  the  comparisons 
of  the  measured  exit  velocities  with  the  predictions  from  the  Model  with  Droplet 
Heat  Transfer  in  Figs.  11.21,  11.22,  and  11.23  are  essentially  a  comparison  of  the 
measured  and  predicted  efficiencies  as  well. 
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11.3  Optimum  Nozzle  Designs 

Based  on  the  parametric  analyses  presented  in  Section  II. 2,  criteria  were 
developed  for  optimum  nozzle  designs.  Desirable  liquid/gas  combinations  were 
investigated  in  terms  of  three  criteria:  maximum  allowable  operating  temperatures 
of  the  liquid,  degradation  characteristics  of  the  liquid,  and  desirable  thermophysical 
properties.  Maximum  efficiencies  for  two-phase  nozzles  were  estimated  based  on 
the  parametric  analyses  of  Section  II. 2.  Finally,  a  simplified  design  procedure 
was  developed  for  the  quick,  first-cut  sizing  by  hand  calculation  of  two-phase 
nozzle  geometry. 

II. 3.1  Criteria  for  Optimum  Nozzle  Designs 

The  parametric  analyses  presented  in  Section  11.2  are  the  bases  for  developing 
the  criteria  for  optimum  nozzle  designs.  First,  considering  the  selection  of 
nozzle  geometry,  nozzle  efficiency  increases  as  nozzle  length  is  increased,  as 
seen  from  Fig.  11.16.  The  increased  length  provides  more  residence  time  for  the 
exchange  of  momentum  and  heat  between  the  two  phases,  and  the  nozzle  efficiency  is 
therefore  higher.  Hence,  within  the  geometric  constraints  of  the  given  application, 
longer  nozzles  result  in  higher  nozzle  efficiency. 

Second,  the  other  governing  flow  parameters  of  the  nozzle  design  are  generally 
set  by  overall  system  considerations.  The  effect  of  these  parameters  on  nozzle 
performance  can  be  studied  for  parametric  trends  that  can  be  used  to  influence  the 
selection  of  the  overall  system  operating  conditions  at  the  design  point.  For 
example,  the  effect  of  increasing  the  nozzle  pressure  ratio,  seen  in  Fig.  II. 6 
(obtained  by  reducing  the  exit  pressure),  is  to  reduce  the  design  point  nozzle 
efficiency.  The  effect  of  nozzle  inlet  pressure  (for  a  given  exit  pressure)  was 
seen  in  Fig.  II. 4,  wherein  lower  nozzle  inlet  pressure  was  seen  to  reduce  the 
nozzle  efficiency  and  can  be  significant  in  the  selection  of  overall  system  operating 
conditions.  In  general,  higher  pressure  ratio  can  be  obtained  by:  increasing  the 
nozzle  inlet  pressure  (this  increase  has  side  effects  on  the  sizing  of  the  pump); 
and/or  by  decreasing  the  exit  pressure  (this  decrease  has  side  effects  both  on  the 
sizing  of  the  turbine  with  the  concomitant  lower  turbine  inlet  pressure  and  on  the 
design  of  the  condenser  where  the  steam  is  condensed  before  being  repressurized  as 
a  liquid). 

The  effect  of  inlet  temperature  on  nozzle  performance  predictions  was  seen 
in  Table  II. 5,  wherein  reducing  the  inlet  temperature  was  seen  to  reduce  the  nozzle 
efficiency  (because  the  droplets  are  larger).  The  effect  of  loading  ratio  was 
seen  in  Fig.  II. 2,  wherein  the  shape  of  the  curve  results  in  a  shallow  minimum 
(instead  of  a  maximum  of  any  kind).  Hence,  loading  ratio  would  probably  be  a 
secondary  factor  from  the  viewpoint  of  nozzle  efficiency  when  selecting  overall 
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system  design  conditions,  unless  a  loading  ratio  lower  than  about  five  was  being 
contemplated  (which  is  the  approximate  value  where  the  efficiency  begins  to  change 
quickly).  Thus,  the  nozzle  performance  trends  that  are  caused  by  the  effects  of 
pressure,  temperature,  and  loading  ratio  can  be  inspected  and  used  to  assist  in 
the  selection  of  the  value  of  these  parameters  for  the  overall  system  conditions. 

The  nozzle  performance,  however,  is  only  one  factor  in  the  overall  design-point 
selection,  because  all  the  components  in  the  system  must  be  considered  collectively 
to  select  the  overall  system  performance. 

II. 3. 2  Liquid/Gas  Combinations 

Suitable  liquid/gas  combinations  were  examined  in  terms  of  three  criteria: 
maximum  allowable  operating  temperatures  of  the  liquid,  degradation  characteristics 
of  the  liquid,  and  desirable  thermophysical  properties.  The  first  two  of  these 
three  criteria  themselves  rule  out  many  possible  combinations,  but  other  important 
criteria  could  include  cost,  availability,  and  toxicity,  depending  upon  the 
application.  High-temperature  organic  heat  transfer  fluids  with  a  low  vapor 
pressure  were  surveyed  as  candidate  liquids  by  examining  reports  in  the  literature 
and  manufacturers'  data  brochures  (which  sometimes  compare  their  own  fluids  with 
those  of  their  competitors).  One  problem  with  some  of  the  fluids  is  that  the 
thermophysical  properties  have  been  measured  only  at  temperatures  near  room 
temperature . 

The  maximum  allowable  operating  temperature  is  an  important  consideration  and 
may  represent  a  limitation  on  the  selection  of  design  conditions.  Table  II. 8 
summarizes  this  maximum  temperature  for  a  number  of  organic  heat  transfer  fluids 
of  possible  interest.  The  degradation  data  in  this  table  were  assembled  from 
research  published  by  Dow  Chemical  (Refs.  10  and  11)  and  show  their  Dowtherm  A 
(DTA)  to  have  relatively  low  degradation  characteristics  at  1060  R;  for  DTA,  their 
extrapolated  measurements  indicate  a  degradation  rate  of  0.0035  weight  percent  per 
week  at  1060  R,  or  approximately  0.2  weight  percent  of  the  DTA  will  degrade  in  one 
year.  This  degradation  is  essentially  a  therma I ly-i nduced  breakdown,  in  a  clean 
closed  system,  of  the  petroleum  oils  into  primarily  volatile  materials.  Furthermore, 
these  degradation  tests  were  conducted  in  a  clean  system,  with  no  oxygen  or  water 
present.  In  a  telephone  conversation  with  Monsanto  (Ref.  12),  their  representative 
indicated  that  Therminol  66  becomes  highly  oxidized  in  the  presence  of  air  or 
water  at  temperatures  higher  than  about  610  R,  and  they  would  not  recommend  the 
use  of  Therminol  66  in  a  two-phase  system  with  steam  as  the  gas  phase.  This 
breakdown  problem  exacerbated  by  the  presence  of  water  could  conceivably  exist  for 
the  other  heat  transfer  oils,  because  "contamination"  by  impurities  is  known  to 
promote  degradation. 

Nevertheless,  based  on  the  information  in  Table  11.8,  DTA  has  a  high  maximum 
operating  temperature  and  relatively  favorable  degradation  characteristics,  pending 
experimental  evidence  for  mixtures  of  DTA  and  steam  at  temperatures  of  interest. 
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Sufficiently  complete  sets  of  thermophysical  data  were  available  only  for 
DTA  and  Dowtherm  G  (DTG),  as  surface  tension  data  (needed  for  estimation  of 
droplet  size)  is  frequently  unavailable,  except  sometimes  for  a  value  at  room 
temperature.  And  the  use  of  Therminol  66  was  definitely  ruled  out  on  the  basis 
of  discussion  with  the  manufacturer's  representative. 

Figure  11.24  compares  the  nozzle  performance  characteristics  for  DTA  and 
DTG.  Basically,  the  main  differences  in  the  thermophysical  properties  at  1060  R 
are  that,  for  DTG,  the  density  is  about  ten  percent  higher,  the  specific  heat  is 
about  seven  percent  lower,  and  the  extrapolated  surface  tension  is  about  42 
percent  higher  than  the  corresponding  properties  of  DTA.  For  loading  ratios  in  the 
range  of  5  to  40,  the  use  of  DTG  instead  of  DTA  would  lead  to  a  one  percentage 
point  drop  in  nozzle  efficiency.  The  exit  temperatures  are  slightly  lower  because 
of  the  lower  liquid  specific  heat,  while  the  exit  liquid  velocities  are  about  one 
percent  lower  and  the  exit  gas  velocities  are  about  one  percent  higher  primarily 
because  of  the  higher  surface  tension  of  DTG.  From  the  viewpoint  of  nozzle 
efficiency,  DTA  is  therefore  better  than  DTG  as  the  liquid-phase  fluid. 

In  the  absence  of  actual  thermophysical  data  for  other  liquids,  the  parametric 
effects  of  the  thermophysical  properties  on  nozzle  performance  can  be  examined  for 
clues  to  be  used  in  searching  for  liquid/gas  combinations  with  favorable  property 
characteristics.  Of  course,  the  effects  of  these  properties  on  nozzle  performance 
may  be  quite  different  from  their  effects  on  overall  system  performance,  and  the 
selected  liquid/gas  combination  must  obviously  have  favorable  overall  system 
performance  characteristics  rather  than  just  favorable  nozzle  performance  character¬ 
istics.  As  seen  earlier.  Fig.  11.13  shows  the  effect  of  liquid  surface  tension  on 
nozzle  performance  predictions.  With  the  other  liquid  properties  being  those  of 
DTA,  a  reduction  in  surface  tension  results  in  an  increase  in  nozzle  efficiency. 
Similarly,  a  decrease  in  liquid  density  as  seen  earlier  in  Fig.  11.14,  or  a 
decrease  in  liquid  specific  heat  as  seen  earlier  in  Fig.  11.12,  will  tend  to 
increase  the  nozzle  efficiency.  As  for  the  effect  of  the  gas-phase  properties,  an 
increase  in  gas  density  as  seen  earlier  in  Fig.  II. 11  will  also  tend  to  increase 
the  nozzle  efficiency,  but  an  increase  in  gas  specific  heat  as  seen  in  Table  II. 6 
will  result  in  only  a  slight  increase  in  nozzle  efficiency.  Thus,  the  trends  in 
nozzle  efficiency  caused  by  the  thermophysical  properties  of  the  liquid  and  of  the 
gas  can  be  used  as  preliminary  guidance  in  examining  candidate  liquids  and  gase6. 

11.3.3  Maximum  Nozzle  Efficiencies 

The  effects  of  all  the  governing  flow  parameters  on  nozzle  efficiency  have 
been  investigated  and  the  results  presented  above  in  Section  II. 2.  Figure  II. 6  is 
a  particularly  relevant  figure  to  examine  for  the  level  of  maximum  nozzle  efficiencies 
that  can  be  expected,  because  those  operating  conditions  are  possibly  typical  of 
those  expected  for  a  two-phase  propulsion  system  with  a  two-phase  turbine.  In 
this  figure,  the  nozzle  inlet  temperatures  are  960  R,  and  the  inlet  pressure  is 
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600  psia.  A  loading  ratio  in  the  range  of  20  to  40  is  quite  reasonable.  Then, 
for  a  nozzle  exit  pressure  of  15  psia,  the  pressure  ratio  is  40,  and  the  efficiency 
of  a  six-inch-long  nozzle  with  a  6team/DTA  mixture  at  a  loading  ratio  of  40  is 
therefore  74  percent.  For  a  steam/DTG  mixture  with  these  operating  conditions, 
the  design  efficiency  would  be  about  one  percentage  point  lower,  based  on  the 
relative  results  of  Fig.  11.24;  a  nozzle  efficiency  of  73  percent  is  therefore 
predicted  with  DTG  as  the  liquid  phase. 


II. 3. 4  Simplified  Procedure  for  Nozzle  Sizing 

Although  the  computerized  model  developed  during  Phase  1  (Ref.  1)  is  relatively 
straightforward  to  use  because  it  is  oriented  toward  engineering  results,  the  code 
still  requires  a  computer  to  obtain  the  results.  As  such,  the  process  to  obtain 
results  is  still  more  complicated  than  a  procedure  for  the  quick  estimation,  by 
hand,  of  first-cut  sizing  numbers  for  the  nozzle,  such  as  the  diameters  at  the 
inlet,  throat,  and  exit.  Hence,  a  simplified  sizing  procedure  was  developed  to 
obtain  quickly  the  approximate  values  of  these  diameters.  This  procedure  is  based 
on  the  isothermal  theory  developed  by  Rudinger  (Ref.  13)  because  this  is  the  most 
complex  theory  that  can  still  be  solved  with  just  a  few  hand  calculations. 

Est_imat_ion_o_f  NozzJ_e_G£ome£r£ 

Five  geometric  parameters  can  be  used  to  approximate  the  nozzle  shape. 

These  parameters  are:  nozzle  length  and  throat  location,  and  the  nozzle  diameters 
(or  flow  areas  if  the  nozzle  is  not  axisymmetric)  at  the  inlet,  throat,  and  exit 
locations.  Based  on  the  parametric  analysis  of  Section  II. 2,  in  order  to  maximize 
nozzle  efficiency,  the  nozzle  length  should  be  as  long  as  possible  within  the 
geometric  constraints  of  the  proposed  application.  Then  the  nozzle  throat  is 
placed  at  the  half-way  point  of  the  nozzle  length,  with  the  throat  pressure 
assigned  a  value  half-way  between  the  inlet  and  exit  pressures.  The  selection  of 
the  throat  location  is  also  based  on  Che  parametric  analysis  of  Section  11.2, 
where  the  flow  area  in  the  vicinity  of  the  throat  was  seen  to  have  a  shallow  shape 
as  a  function  of  nozzle  length  and  where  the  location  of  the  throat  generally 
occurred  near  the  mid-point  of  the  nozzle.  Hence,  for  the  purpose  of  the  first-cut 
sizing  calculations,  these  assumptions  pin  down  two  of  the  geometric  parameters  of 
the  nozzle  design:  the  nozzle  throat  location  and  the  nozzle  length. 

Estimation  of  Inlet  Diameter 


The  inlet  diameter  (or  flow  area)  is  estimated  using  the  continuity  equation 
and  assuming  that  the  total  flow  area  is  the  sum  of  the  flow  areas  required  for 
each  of  the  two  phases: 
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The  quantities  needed  to  estimate  the  nozzle  inlet  flow  area  are  all  known  from 
the  assumed  inlet  conditions,  and  the  perfect  gas  law  is  used  to  calculate  the  gas 
density.  For  the  propulsion  application,  the  velocity  ratio  (V^/Vg)  is  generally 
in  the  range  of  0.6  to  0.8,  so  a  value  of  0.7  is  generally  a  good  choice  to  be  used 
for  the  first-cut  sizing  calculations.  Hence,  the  nozzle  inlet  diameter  (d^n)  can 
bt  estimated  from  the  inlet  flow  area: 
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Estimation  of  Throat  Diameter 


The  nozzle  throat  diameter  is  estimated  by  using  the  continuity  equation  in 
conjunction  with  the  velocity  equation  which  results  from  the  momentum  balance  in 
the  isothermal  theory  developec  by  Rudinger  (Ref.  13).  The  gas  velocity  at  the 
nozzle  throat  is  estimated  with  the  expression: 


(11.13) 


where  r  is  the  loading  ratio,  K  is  the  velocity  ratio  (Vj/Vg),  pthr  *8  the 
pressure  at  the  throat  and  is  assumed  to  be  half  way  between  the  inlet  and  exit 
pressures,  R  is  the  gas  constant,  and  e^n  is  the  inlet  value  of  the  liquid 
volume  fraction: 
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The  velocity  ratio  is  assumed  constant  at  the  inlet  value,  because  this  assumption 
was  made  by  Rudinger  in  order  to  obtain  a  closed-form  analytical  solution.  The 
throat  liquid  velocity  is  calculated  from  the  velocity  ratio:  Vt  thr  *  K  Vg 
Then  by  continuity,  the  throat  diameter  (dt^r)  is  estimated  with  the  expression: 
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where  the  liquid  volume  fraction  at  the  throat  is  estimated  with: 
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The  gas  density  at  the  throat, P^  tjlr  is  calculated  from  the  perfect  gas  law 
with  the  assumed  throat  pressure  and  from  the  inlet  value  of  the  gas  temperature 
because  the  flow  is  assumed  to  be  isothermal. 


Estimation  of  Exit  Diameter 


A  similar  approach  is  used  to  estimate  the  nozzle  exit  diameter.  The  gas 
velocity  at  the  nozzle  exit  is  estimated  with  the  momentum-derived  expression: 
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Then  the  nozzle  exit  diameter  is  estimated  with  the  continuity  expression: 
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where  the  liquid  volume  fraction  is  estimated  with: 
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The  gas  temperature  and  both  the  loading  ratio  and  the  velocity  ratio  are  constant 
throughout  the  nozzle,  as  per  the  original  assumptions  for  the  Isothermal  Theory. 
Thus,  the  nozzle  diameters  at  the  three  locations  of  especial  interest  (inlet, 
throat,  and  exit)  can  be  estimated  with  Eqs.  11.12,  11.15,  and  11.18,  respectively. 
Table  11.9  summarizes  the  equations  for  the  Simplified  Nozzle  Sizing  Procedure. 
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£omparj_s£n_ojf  £a_l_C£l_at_io_ns 

Calculations  were  made  to  obtain  an  indication  of  the  differences  that 
might  be  expected  between  the  diameters  from  the  Simplified  Nozzle  Sizing 
Procedure  and  the  diameters  obtained  from  the  computerized  Model  with  Droplet 
Heat  Transfer.  Figure  11.25  compares  the  results  for  nozzle  geometry  from  the 
two  models.  At  these  operating  conditions,  the  results  for  nozzle  diameter 
from  the  simplified  sizing  procedure  at  a  velocity  ratio  of  either  0.7  or  0.8 
are  within  eight  percent  of  the  more  exact  results  from  the  computerized  model. 
The  results  for  nozzle  efficiency,  obtained  by  inserting  the  velocity  results 
into  the  definition  of  nozzle  efficiency  (Eq.  11.9),  also  show  reasonable 
agreement,  except  at  the  lower  end  of  the  loading  ratio  where  the  assumption  of 
isothermal  conditions  throughout  the  nozzle  is  least  defendable.  The  drop  in 
gas  temperature  increases  as  the  loading  ratio  is  reduced  (as  seen  earlier  in 
Fig.  11.2)  in  the  actual  situation,  whereas  increased  amounts  of  heat  addition 
would  be  required  if  the  isothermal  assumption  were  to  be  maintained  in  the 
Isothermal  Model.  Furthermore,  as  seen  in  Fig.  11.26,  the  exit  velocities  of 
the  liquid  that  are  calculated  with  the  Simplified  Sizing  Procedure  are  lower 
than  those  calculated  with  the  computerized  Model  with  Droplet  Heat  Transfer, 
and  the  liquid  phase  carries  the  bulk  of  the  kinetic  energy  for  loading  ratios 
greater  than  about  three  (as  seen  earlier  in  Fig.  11.17).  Thus,  a  lower 
calculated  liquid  exit  velocity  explains  the  lower  calculated  values  of  nozzle 
efficiency  from  the  Simplified  Nozzle  Sizing  Procedure.  Further,  because  the 
exit  value  of  the  velocity  ratio  is  actually  closer  to  0.8  for  the  operating 
conditions  of  Fig.  11.26,  the  nozzle  efficiencies  in  Fig.  11.25  that  were 
calculated  from  the  Simplified  Sizing  Procedure  at  a  velocity  ratio  of  0.8  are 
much  closer  to  the  more  exact  results  for  nozzle  efficiencies  from  the  Model 
with  Droplet  Heat  Transfer. 

Figure  11.27  compares  the  results  for  nozzle  efficiency  and  diameters  from 
the  two  models  at  a  different  set  of  operating  conditions.  Here,  at  a  velocity 
ratio  of  0.7,  the  results  from  the  Simplified  Sizing  Procedure  for  the  throat 
diameter  are  about  four  percent  lower  than  the  results  from  the  Droplet  Heat 
Transfer  Model,  and  the  nozzle  efficiency  results  agree  within  two  percent 
except  for  loading  ratios  lower  than  about  five,  while  the  results  for  the  exit 
diameter  are  about  18  percent  higher  than  those  from  the  Droplet  Heat  Transfer 
Model.  Thus,  the  comparison  of  results  in  Figs.  11.25  and  11.27  are  an  indica¬ 
tion  of  the  differences  that  can  be  expected  between  the  quick  first-cut 
calculations  of  the  Simplified  Nozzle  Sizing  Procedure  and  the  more  exact 
results  of  the  Model  with  Droplet  Heat  Transfer. 
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II .4  Assessment  of  Work  Extraction  Potential 

In  this  section,  performance  criteria  of  two-phase  nozzles  applicable  to 
work  extraction  systems  are  identified,  and  the  attributes  of  a  nozzle  suitable 
for  inclusion  in  a  system  that  converts  kinetic  energy  into  work  are  assessed 
in  terms  of  the  off-design  performance  characteristics  of  a  two-phase  nozzle. 

These  off-design  characteristics  are  also  of  interest  because  in  real  applications, 
the  nozzle  will  most  certainly  operate  at  conditions  other  than  the  design-point 
conditions  for  which  it  was  designed. 

11.4.1  Performance  Criteria  for  Work  Extraction 

One  of  the  two-phase  nozzle  performance  criteria  applicable  to  work 
extraction  systems  is  the  nozzle  efficiency,  which  is  defined  as  the  increase 
in  kinetic  energy  of  the  two-phase  mixture  relative  to  the  ideal  change  in 
enthalpy  if  the  two-phase  mixture  were  to  be  expanded  isentropical ly  as  a 
homogeneous  mixture  to  the  same  exit  pressure.  The  kinetic  energy  of  the 
two-phase  mixture  at  the  nozzle  exit  is  in  turn  converted  into  work,  as  seen  in 
the  schematic  flow  diagram  of  Fig.  11.28.  Thus,  the  maximum  nozzle  efficiency 
corresponds  to  the  maximum  level  of  conversion  to  kinetic  energy  that  is 
possible  by  the  nozzle  at  the  assumed  operating  conditions.  At  the  nozzle 
inlet,  the  heated  liquid  (DTA ,  for  example)  is  mixed  with  the  condensate  stream 
(water,  for  example)  to  produce  the  gas  phase  by  flashing  the  condensate  stream 
if  the  enthalpy  of  the  heated  liquid  is  sufficiently  high.  Then  the  mixture  is 
expanded  and  the  enthalpy  of  the  mixture  is  converted  into  kinetic  energy  for 
subsequent  conversion  to  shaft  power  in  the  two-phase  turbine.  As  seen  in  Fig. 
11.17,  most  of  the  kinetic  energy  is  carried  by  the  liquid  stream  over  the  range 
of  loading  ratios  which  are  of  interest  for  propulsion  applications.  Hence  in 
order  to  extract  most  of  the  mixture's  energy,  the  liquid-phase  kinetic  energy 
must  obviously  be  extracted. 

A  second  nozzle  performance  criterion  applicable  to  work  extraction 
systems  is  the  mixture  velocity  at  the  nozzle  exit.  For  the  shaft  speed  of  the 
two-phase  turbine  to  be  advantageously  lower  than  that  of  an  all-gas  turbine, 
the  nozzle  exit  velocity  must  be  reasonably  low.  Otherwise,  the  low-speed 
advantage  of  a  two-phase  turbine  will  not  be  gained,  and  a  large  gearbox  would 
still  be  required  to  reduce  the  turbine  shaft  speed  down  to  the  thruster  shaft 
speed.  As  presented  in  Section  II. 2,  high  nozzle  exit  velocities  occur  when 
the  loading  ratio  or  the  nozzle  exit  pressure  is  very  low.  Hence,  very  low 
levels  of  loading  ratio  or  nozzle  exit  pressure  should  be  avoided  if  the 
advantage  of  low  turbine  shaft  speed  is  to  be  realized  in  the  work  extraction 
system. 
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The  most  significant  factors  for  the  work  extraction  system  are  the 
overall  system  efficiency  and  design,  not  just  the  nozzle  efficiency  which 
itself  might  be  favorably  enhanced  by  factor(s)  that  could  ultimately  result  in 
a  lower  system  efficiency,  or  a  less  attractive  system  design.  As  discussed 
above,  the  two-phase  nozzle  is  only  one  of  several  components  in  a  work  extrac¬ 
tion  system.  For  example,  a  lower  nozzle  pressure  ratio  (obtained  by  increasing 
the  nozzle  exit  pressure)  would  result  in  a  higher  nozzle  efficiency,  but  the 
effect  of  this  lower  nozzle  pressure  ratio  on  the  overall  system  efficiency 
must  also  be  examined.  Generally,  as  the  nozzle  exit  pressure  increases,  the 
mixture  velocity  at  the  nozzle  exit  will  decrease,  the  turbine  speed  will  be 
lower  and  lower  cycle  efficiency  may  result  (Ref.  14).  Hence,  the  favorable 
reduction  in  gearbox  size  that  is  derived  from  the  lower  nozzle  exit  velocity 
must  be  balanced  against  any  potential  reduction  in  cycle  efficiency  of  the 
work  extraction  system,  and  this  trade-off  optimization  may  depend  in  part  on 
the  particular  application. 

II. 4. 2  Nozzle  Attributes  for  Energy  Conversion  Systems 

In  order  to  assess  the  attributes  of  a  nozzle  suitable  for  inclusion  in  a 
work  extraction  system,  a  simplified  calculation  procedure  was  developed  to 
estimate  the  off-design  performance  characteristics  of  two-phase  nozzles.  The 
off-design  performance  characteristics  are  important  in  assessing  the  attractive¬ 
ness  of  the  two-phase  nozzle  over  the  entire  range  of  expected  operating 
conditions  because  the  overall  system  will  undoubtedly  have  to  operate  at 
conditions  quite  different  from  the  design-point  conditions.  As  in  the  develop¬ 
ment  of  the  Simplified  Nozzle  Design  Procedure,  the  isothermal  theory  of 
Rudinger  (Ref.  13)  was  again  adapted  for  the  off-design  calculations  because 
this  is  the  most  complex  theory  that  is  still  tractable  in  terms  of  doing 
calculations  reasonably  quickly  by  hand. 

Continuity  equations  and  the  momentum  equation  are  used  to  estimate  the 
off-design  nozzle  performance. 

C°n^  inu^J^_R  e^i  reme  n  t  s 

At  a  given  set  of  operating  conditions,  the  continuity  equation  for  the 
gas  flow  can  be  written  for  flow  conditions  anywhere  in  the  nozzle: 


\  ■  f-'ln)  »„ln  Yin  ‘in  "  (l“)  P,», 


(11.20) 


where  c  is  the  local  liquid  volume  fraction  and  is  calculated  with  the  expression: 

1 _ 

(11.21) 


re. 
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Again  K  is  the  velocity  ratio  (Vt/V^)  and  is  assumed  to  be  constant  throughout 
the  nozzle  at  a  given  set  of  operating  conditions  (as  per  the  original  assumption 
that  Rudinger  made  in  order  to  obtain  the  version  of  the  momentum  equation  that 
will  be  used  later).  Hence  from  Eq.  11.20,  the  gas  velocity  at  any  location  in 
the  nozzle,  relative  to  the  inlet  gas  velocity,  is  then  known  as: 


PR» in  Ain 
Pg  A 


(11.22) 


where  the  perfect  gas  law  is  used  to  convert  density  into  pressure,  for  the 
assumption  of  isothermal  conditions  in  the  nozzle  at  one  set  of  operating 
conditions . 

A  second  group  of  continuity  equations  is  used  to  describe  the  gas  flow 
rates  at  different  possible  operating  conditions  at  the  nozzle  inlet,  where  the 
superscript  (*)  denotes  the  design-point  (or  reference)  operating  conditions: 


mg  -  ^  Cinj  pg,in  Vg,in  Ain 

mg  "  (1-£in)  PB«in  Vg.i« 


(11.23) 

(11.24) 


These  equations  can  be  combined  to  obtain  the  inlet  gas  velocity  at  any  operating 
conditions,  relative  to  the  inlet  gas  velocity  at  the  design-point  conditions: 


V 


g.in 


V* 

g ,  in  nig 


(11.25) 


If  the  assumption  is  made  that  the  inlet  conditions  of  temperature  and  pressure 

are  the  same  at  both  design  and  off-design,  then  p  •  *  p  •  ,  and  Eq.  11.25 

.  ,  .  , .  g  > 1 n  S » ln 

simplifies  to: 


(II. 25a) 


Further,  if  the  special  operating  conditions  of  the  same  gas  flow  also  obtain 
at  both  design  and  off-design  conditions,  then  Eq.  II. 25a  further  reduces 
to: 
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V 


g.in 


(II. 25b) 


Equations 

Now , 
Eq.  11.25 


II. 25a  and  II. 25b  obviously  represent  special  circumstances. 


to  calculate  the  gas  velocity  throughout  the  nozzle  at  off-design, 
is  combined  with  Eq.  11.22  to  obtain: 


(11.26) 


which  is  valid  for  all  combinations  of  inlet  conditions  and  for  all  variations 
in  gas  flow  rates.  However,  if  the  assumption  is  made  that  the  inlet  conditions 
of  temperature  and  pressure  are  the  same  for  all  operating  conditions,  then 
4, in  =  eg, in-  And  b*  the  Perfect  gas  law,  pg  in/pg  -  Pin/P  -  PR, 
where  PR  is  the  nozzle  pressure  ratio  at  that  location  in  the  nozzle  where 
the  flow  area  is  A.  Then,  using  Eq .  11.21  for  the  definition  of  e,  Eq.  11.26 
can  be  rewritten  as: 


tr* 


PR 


in 


g.in 


* 


m 


g 


1  + 


r  pg,m 

(PR)Kp5 


1  + 


r*P*  , 

g.in 

K*P„ 


(II. 26a) 


which  is  valid  for  off-design  operation  at  unchanged  values  of  the  inlet 
temperature  and  pressure,  but  at  various  gas  flow  rates  and  loading  ratios. 
Equation  II. 26a  expresses  the  off-design  gas  velocity  (V^),  *t  the  location 
where  the  nozzle  flow  area  is  A,  in  terms  of  the  gas  velocity  (V_  £n)  at 
the  nozzle  inlet  at  the  design-point  conditions.  If  the  nozzle  flow  area  (A) 
is  selected  to  be  the  exit  flow  area,  then  the  gas  velocity  (V^)  is  the 
nozzle-exit  gas  velocity. 


But  still  another  expression  is  required  in  order  to  obtain  the  nozzle 
exit  pressure.  After  all,  when  the  nozzle  was  designed,  the  inlet  conditions 
and  the  pressure  ratio  were  selected.  Now,  at  off-design,  the  inlet  conditions 
are  again  selected,  but  the  exit  area  is  fixed,  and  the  exit  pressure  (hence 
the  pressure  ratio)  must  then  be  calculated  to  be  consistent  with  the  requirements 
of  the  continuity  and  momentum  equations.  The  continuity  requirements  have  been 
developed  above,  and  now  the  equation  to  satisfy  the  momentum  requirements  will 
be  developed. 


Momentum  Requirements 


The  expression  for  gas  velocity  at  the  location  in  the  nozzle  where  the 

pressure  is  P  was  obtained  for  the  isothermal  theory  by  Rudinger  (Ref.  13), 
from  momentum  considerations.  This  expression  is: 
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v:  -  V 


RT 


g.in 


g.in 


1+rK 


(11.27) 


where  R  is  the  perfect  gas  constant.  In  addition  to  the  isothermal  assumption 
discussed  above,  the  assumption  of  constant  velocity  ratio  (K)  throughout  the 
nozzle  at  a  given  set  of  operating  conditions  was  made  by  Rudinger  to  obtain 
Eq.  11.27. 


From  the  continuity  relationship  of  Eq.  11.22,  for  a  given  set  of  operating 
conditions,  the  gas  velocity  at  any  location  is  related  to  the  inlet  gas 
velocity,  and  Eq.  11.22  can  then  be  substituted  into  Eq.  11.27  to  obtain  the 
result : 


2RT  , 

8,  in 

1+rK 

in  | 

M  *  (i 

»..)(■  -  U 

1 

rpin 

Ain  ^in) 
A  1—  c  j 

I2  - 

(11.28) 


Ca_l_cu^at^cm  Procedure^ 

Now  that  the  continuity  and  momentum  relations  in  their  adapted  form  are 
available  for  the  simplified  off-design  performance  predictions,  the  calculation 
procedure  can  be  discussed. 

The  off-design  gas  inlet  velocity  is  calculated  from  Eq.  11.25,  which  is 
an  expression  based  solely  on  continuity.  Here,  the  off-design  conditions  (of 
inlet  gas  velocity,  inlet  liquid  volume  fraction,  velocity  ratio,  loading 
ratio,  gas  mass  flow  rate,  etc.)  relative  to  the  design-point  values  are  used. 
The  off-design  velocity  ratio  must  be  assumed,  and  one  possible  estimation 
method  (that  was  used  here)  is  to  use  the  values  from  the  design-point  calcula¬ 
tions  with  the  Model  with  Droplet  Heat  Transfer  at  other  values  of  the  loading 
ratio  (for  example,  the  series  of  design-point  values  of  Vj/Vg  in  Fig.  II. 2 
at  the  inlet  conditions  of  temperature  and  pressure  of  interest). 

Then,  because  Eq.  11.28  cannot  be  solved  explicitly  for  the  exit  pressure 
ratio  (P^n/P)»  a  trial-and-error  solution  is  required  to  obtain  the  value  of 
Pin/P  that  yields  the  value  of  Vg  £n  that  agrees  with  the  value  from  Eq. 

11.25.  (The  pressure  term  (P)  occurs  in  four  places,  including  the  e  term,  in 
the  right-hand  side  of  Eq.  11.28,  and  an  explicit  solution  for  P  does  not  seem 
possible.)  Then,  knowing  the  values  of  the  pressure  ratio  (Pjn/P)  and  the 
inlet  gas  velocity  (V^  £n)  that  satisfy  Eqs.  11.25  and  11.28  at  the  off-design 
conditions,  the  exit  gas  velocity  can  be  calculated  with  Eq.  11.27. 

With  the  off-design  velocities  at  the  nozzle  inlet  and  exit  now  known,  the 
conversion  of  enthalpy  into  kinetic  energy  by  the  nozzle  can  then  be  estimated. 
The  change  in  kinetic  energy  of  the  two-phase  mixture  as  it  passes  through  the 
nozzle  is  given  by  the  expression: 
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where  m  is  the  total  mass  flow  rate  (m  •  m^  +  m^)  and  7  is  the  mass-weighted 
velocity:  V  *  (V^  +  r  V ^ ) / ( 1  +  r).  Recalling  that  the  velocity  ratio  in  the 
isothermal  theory  is  assumed  to  be  constant  throughout  the  nozzle  at  a  given 
operating  point,  then  Eq.  11.29  can  be  rewritten  as: 


AKE 


1  • 

Img 


Or+rK)2  Ll  _ 
1+r  \  8 


(II. 29a) 


Then  the  kinetic  energy  conversion  at  off-design  relative  to  the  design-point 
conversion  is: 


AKE 

AKE* 


(11.30) 


And  if  the  gas  flow  rate  is  the  same  at  both  design  and  off-design,  then  the 
flow  rate  factor  in  Eq.  11.30  is  unity. 


Calculated  Results 


Results  calculated  from  the  simplified  procedure  developed  above  were  used 
to  investigate  the  effect  of  off-design  operation  on  the  nozzle  performance 
characteristics.  For  the  purpose  of  examining  off-design  trends,  the  inlet 
conditions  of  temperature  and  pressure  were  kept  constant  at  their  design-point 
values.  In  actual  operation,  the  nozzle  inlet  temperature  and  pressure  could 
vary  dramatically  at  the  off-design  operating  conditions,  depending  on  the 
operation  of  the  other  components  in  the  work  extraction  system  (Fig.  11.28)  at 
off-design,  but  these  off-design  inlet  conditions  were  arbitrarily  chosen  because 
no  actual  schedule  was  available.  Figure  11.29  shows  results  for  the  off-design 
performance  of  the  two-phase  nozzle,  assuming  that  the  gas  flow  rate  is  kept  constant 
at  the  design-point  value  while  the  liquid  flow  rate  (and  hence  the  loading  ratio) 
is  varied.  Other  schedules  of  gas  flow  rate  could  easily  be  chosen,  but  this 
schedule  was  chosen  for  this  set  of  calculations.  The  velocity  ratio  selected 
was  based  on  the  design-point  values  over  the  range  of  loading  ratio  calculated 
previously  with  the  Model  with  Droplet  Heat  Transfer.  As  the  loading  ratio  is 
reduced  in  Fig.  11.29,  the  inlet  gas  velocity  is  seen  to  decrease,  as  required 
by  continuity  considerations  for  the  same  inlet  conditions  (this  is  due  to  the 
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fact  that  less  mass  is  going  through  the  nozzle  inlet),  while  the  exit  gas 
velocity  increases  and  the  nozzle  pressure  ratio  increases  (actually  the  exit 
pressure  decreases,  because  the  nozzle  inlet  pressure  is  held  constant  in  these 
calculations).  With  a  reduced  loading  ratio  but  at  a  fixed  gas  flow  rate,  the 
gas  will  expand  more  because  less  liquid  mass  is  present  to  retard  the  gas 
expansion.  Generally,  as  seen  earlier  in  the  design-point  calculations  (for 
example.  Fig.  11.2),  as  the  loading  ratio  is  reduced,  the  gas  velocity  increases 
because  the  gas  expands  more  in  the  presence  of  less  liquid.  Further,  to  allow 
the  increased  gas  expansion,  the  nozzle  exit  pressure  will  be  lower  (hence,  the 
gas  density  will  be  lower,  to  satisfy  continuity  for  the  isothermal  expansion). 
The  relative  kinetic  energy  conversion,  calculated  with  Eq.  11.30,  shows  an 
increase  of  less  than  five  percent  above  the  design-point  value.  Similar 
results  are  shown  in  Fig.  11.30  for  different  inlet  conditions,  and  with  the 
same  assumption  of  the  same  gas  flow  rate  at  design  as  at  off-design  conditions. 
Here,  the  expected  rang*  of  the  velocity  ratio  as  a  function  of  loading  ratio 
was  assumed  to  be  wider  (based  on  previous  design  point  calculations  from  the 
Model  with  Droplet  Heat  Transfer,  at  the  lower  values  of  inlet  temperature  and 
pressure).  This  larger  range  may  explain  why  the  effect  of  loading  ratio  on 
the  relative  kinetic  energy  conversion  is  slightly  larger  while  overall  it  is 
still  a  small  effect  in  that  it  varied  by  less  than  eight  percent  over  a  range 
of  loading  ratio  from  double  to  half  the  design-point  value. 

Another  possible  mode  of  off-design  operation  is  for  the  oading  ratio  to 
be  kept  constant  over  the  entire  range  of  operating  conditions.  In  this  mode, 
both  the  liquid  flow  rate  and  the  gas  flow  rate  are  varied  in  tandem  to  maintain 
a  constant  value  of  the  loading  ratio.  Figure  11.31  shows  the  results  calculated 
for  6uch  a  mode,  wherein  the  inlet  conditions  of  temperature  and  pressure  were 
again  held  constant  at  the  design-point  values,  throughout  the  off-design 
operation.  As  the  gas  flow  rate  decreases  (while  still  at  constant  loading 
ratio),  the  inlet  gas  velocity  decreases  because  less  total  mass  is  flowing 
through  the  nozzle,  while  the  exit  gas  velocity  is  increasing.  The  inlet  gas 
velocity  is  lower  because  of  continuity  (Eq.  11.25);  while  to  satisfy  the 
momentum  equation  in  the  form  of  Eq.  11.28,  the  pressure  ratio  must  be  increased 
(by  reducing  the  exit  pressure  in  this  instance  because  the  inlet  pressure  is 
held  constant)  and  an  increased  pressure  ratio  means  a  higher  exit  gas  velocity, 
Eq.  11.27.  Simplified,  in  terms  of  the  isothermal  theory,  a  lower  total  mass 
flow  through  a  nozzle  that  is  designed  for  a  higher  mass  flow,  will  be  expanded 
further  than  the  higher  mass  flow,  for  these  levels  of  the  loading  ratio. 

The  relative  kinetic  energy  conversion  for  the  operating  mode  of  constant 
loading  ratio  is  also  shown  in  Fig.  11.31  and  is  seen  to  decrease  as  the 
loading  ratio  is  decreased,  because  the  lower  mass  flow  rate  more  than  makes  up 
for  the  higher  exit  gas  velocity.  Actually  the  increase  in  gas  exit  velocity 
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is  relatively  small  as  compared  to  the  increase  seen  in  Fig.  11.30  for  the  mode 
of  operation  where  the  loading  ratio  varies.  This  mode  (with  the  constant  gas 
flow  rate  over  the  off-design  conditions)  has  a  much  flatter  profile  of  kinetic 
energy  conversion  at  off-design  than  for  the  mode  where  the  gas  flow  rate  is 
varied  (and  the  loading  ratio  is  held  constant).  This  difference  could  possibly 
be  useful  when  the  off-design  conditions  for  the  work  extraction  system  are 
defined  because  the  liquid  flow  rates  and  the  gas  flow  rates  could  conceivably 
be  scheduled  independently  to  obtain  the  flatter  profile. 

Thus,  using  the  simplified  calculation  procedure  described  here,  and 
within  the  limitations  and  assumptions  of  this  procedure,  the  off-design 
performance  of  two-phase  nozzles  can  be  estimated.  Consequently,  the  off-design 
flow  rates  might  be  tailored  to  vary  the  shape  of  the  off-design  profile  of 
kinetic  energy  conversion  by  the  nozzle,  for  appropriate  integration  of  the 
nozzle  into  an  energy  conversion  system.  A  valuable  design-procedure  outgrowth 
of  these  studies  is  that,  because  the  exit  gas  velocity  for  the  two-phase 
nozzle  does  decrease  as  the  mass  flow  is  increased,  consideration  should  be 
given  to  selecting  the  design-point  condition  for  the  nozzle  as  a  condition 
where  the  mass  flow  is  relatively  low  in  comparison  to  the  overall  schedule  of 
mass  flows,  so  that  the  nozzle  exit  velocity  will  be  lower  at  the  other  expected 
(off-design)  operating  points  where  the  mass  flow  rates  would  be  higher.  This 
selection  would  then  tend  to  retain  the  advantage  of  low  nozzle  exit  velocities 
throughout  the  entire  expected  range  of  nozzle  operation. 
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TABLE  11.1 

SENSITIVITY  OF  NOZZLE  PERFORMANCE  PREDICTIONS  TO  DROPLET  BREAKUP 


SIX-INCH  NOZZLE  LENGTH 
STEAM  DTA  MIXTURE 
INlET  1060  R  ANL  lOL'L1  ps'd 
EXIT  lOL'DSd 


VALUES  AT  EXIT 


1 

We)cnt 

V9 

V{ 

Vf/Vg 

T9 

T| 

D 

V 

ft/sec 

R 

ft 

1 

6 

2570 

2226 

0  866 

8066 

924.1 

2.23  xl  O'5 

0  930 

1 

3 

2554 

2272 

0.890 

811  8 

914.7 

1.71  X  10'5 

0  942 

1 

1 

2532 

2331 

0  921 

818.3 

902.8 

1.16X10'5 

0  957 

5 

6 

1661 

1381 

0.831 

942  6 

1004.6 

3.19  x  10'5 

0883 

5 

3 

1633 

1406 

0  861 

948  3 

1001.7 

2  .46  x  1 0'5 

0  903 

5 

1 

1598 

1438 

0  890 

954.9 

998.0 

1.68  x  1 0'5 

0  930 

10 

6 

1310 

1064 

0812 

985  7 

1028  2 

4.03x1  O'5 

0  873 

10 

3 

1281 

1082 

0  845 

990.3 

1026  7 

3.11  x  10‘5 

0  895 

10 

1 

1244 

1104 

0  888 

9956 

1024.7 

2  12  x  10'5 

0  923 

40 

6 

837  1 

653  4 

0  781 

1031  3 

1051  0 

7.04  x  10‘5 

0  890 

40 

3 

809  5 

661  0 

0817 

1034  0 

1050  7 

5.42  XlO’5 

0  908 

40 

1 

775  2 

671  0 

0  866 

1037  1 

1050  2 

3.68  x  10*5 

0  933 

80-11-64-3 
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TABLE  11.2 

SENSITIVITY  OF  NOZZLE  PERFORMANCE  PREDICTIONS  TO  DROPLET  DRAG 

SIX  INCH  NOZZLE  LE  NGTH 
STEAM  DTA  MIXTURE 
INLET  T'JoO  R  AND  <000  PS’d 
EXIT  lOL’PS-rf 


VALUES  AT  EXIT 


1 - 

r 

CD 

V9 

Vf 

Vf/Vn 

Tg 

D 

V 

% 

ft/sec 

i  g 

R 

(t 

i 

80 

2579 

2216 

0  859 

812.1 

918  4 

2  05  x  10"  5 

0  930 

i 

100 

2570 

2226 

0.866 

806  7 

924  1 

2  24  xl0'5 

0  930 

i 

120 

2562 

2233 

0  872 

802  0 

929  1 

2.40  x  1 0'5 

0  930 

5 

80 

1671 

1376 

0  823 

949  1 

1003.5 

2.91  x  10'5 

0  880 

5 

100 

1661 

1381 

0  831 

942.6 

1004  6 

3.19  x  10"  5 

0  883 

5 

120 

1653 

1384 

0.837 

936  8 

1005  6 

3  44  X  10-5 

0  885 

10 

80 

1320 

1060 

0.804 

990  0 

1027.8 

3.67x1 0’5 

0  869 

10 

100 

1310 

1064 

0  812 

985.7 

1028  2 

4  03  X  10’ 5 

0.873 

10 

120 

1303 

1066 

0.819 

980  8 

1028  5 

4.35x1 0’5 

0  876 

40 

80 

844  7 

651.7 

0.772 

1034  2 

1051.0 

6  41  x  10-5 

0  886 

40 

100 

837.1 

653  4 

0  781 

1031  3 

1051  0 

7  05  x  10-5 

0  890 

40 

120 

831  2 

654  7 

0.788 

1028  6 

1051.0 

7.62  x  10*5 

0  893 
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TABLE  11.3 

SENSITIVITY  OF  NOZZLE  PERFORMANCE  PREDICTIONS  TO  HEAT  TRANSFER 

SIX-INCH  NOZZLE  LENGTH 
STEAM.'DTA  MIXTURE 
INLET  1  060  R  AND  1000  psia 
EXIT  lOOpsia 


VALUES  AT  EXIT 


r 

HTC 

V9 

Vf/Vg 

T9 

Tf 

T? 

% 

ft/sec 

R 

i 

80 

2564 

2224 

0.867 

799  4 

932.7 

0  927 

i 

100 

2570 

2226 

0.866 

806.6 

924.1 

0.930 

i 

120 

2575 

2227 

0.865 

812.1 

917.5 

0933 

5 

80 

1657 

1378 

0832 

9336 

1006  7 

0  880 

5 

100 

1661 

1381 

0  831 

942  6 

1004  6 

0  883 

5 

120 

1664 

1382 

0  830 

9492 

1003  0 

0886 

10 

80 

1308 

1063 

0.813 

978  1 

1029  0 

0  871 

10 

100 

1310 

1064 

0.812 

985  7 

1023  2 

0  873 

10 

120 

1312 

1065 

0.811 

991.0 

1027.6 

0  875 

40 

80 

836  5 

653  2 

0  781 

1027  0 

1051  1 

0889 

40 

100 

637  1 

653  4 

0781 

1031  3 

1051  0 

0  890 

40 

120 

837  9 

653  7 

0  780 

1034.3 

10509 

0  890 

1 


TABLE  11.4 

RELATIVE  EFFECT  OF  GOVERNING  PARAMETERS  ON  NOZZLE  PERFORMANCE 
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NOZZLE  LENGTH 
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TABLE  11.5 

EFFECT  OF  INLET  TEMPERATURE  ON  NOZZLE  PERFORMANCE  PREDICTIONS 

STtAM  D1A 
Six. INCH  NOZZi-t 

FLCU\  CONDITIONS  P„  =  600  Pb  re>  ,  =-  IOC'  PS- 


R81 -955229-4 


TABLE  11.6 

EFFECT  OF  GAS  SPECIFIC  HEAT  ON  NOZZLE  PERFORMANCE  PREDICTIONS 


GAS  HIA 
Six  iNC'h  noZZlE 

FlOv\  C ON L'  ■  t IONS  fqi,n  =  -  lutC  H  P(r  =  ’000  pb'a 

Pe>  .  =  UK'  ps'd 


LOADING 

RATIO 

r 

THERMAL 
CAPACITY 
RATIO 
(C  CDgl,n  ' 

r  -  fcAl  1  VALUto - - - - - -  ’  ] 

V9 

T9 

Tf 

7? 

tisec 

R 

c 

40 

i  ’  8 

837  e 

653  7 

1033  5 

1051  6 

0.8902 

1 

4; 

8-  A 

65  j  5 

1031  3 

1051  0 

08898 

0  5 

40 

7 «  . 

- 

653  4 

1029  9 

1050  7 

0  8896 

2 

10 

4  ‘ 

13U 

1066 

997  9 

1033  0 

0875 

1 

10 

e  9 

1310 

1064 

985  7 

1028  2 

0  873 

0  5 

10 

1  7  8 

1308 

1062 

977  1 

1025  0 

0.871 

b 


r 

vg)in 

Vf  )in 

n 

VfA/g 

T9 

Tf 

V 

ft/sec 

<t/sec 

R 

10 

100 

70 

1310 

1064 

0.812 

986 

1028 

0.873 

10 

200 

140 

1320 

1073 

0.813 

985 

1028 

0.876 

10 

300 

210 

1334 

1087 

0.815 

985 

1028 

0.880 

20 

100 

70 

1035 

822 

0.794 

1014 

1043 

0.874 

20 

200 

140 

1046 

833 

0.796 

1014 

1043 

0878 

20 

300 

210 

1064 

850 

0.799 

1013 

1043 

0.883 

40 

100 

70 

837 

654 

0.781 

1031 

1051 

0.890 

40 

200 

140 

850 

666 

0.784 

1031 

1051 

0.894 

40 

300 

210 

871 

687 

0.789 

1030 

1051 

0.899 

20 

100 

70 

1035 

822 

0.794 

1014 

1043 

0.874 

20 

100 

80 

1036 

823 

0  794 

1014 

1043 

0.875 

20 

100 

90 

1038 

824 

0.794 

1014 

1043 

0.875 

20 

300 

210 

1064 

850 

0.799 

1013 

1043 

0.883 

20 

300 

240 

1072 

858 

0.801 

1013 

1043 

0  885 

20 

300 

270 

1081 

868 

0.802 

1012 

1043 

0.887 
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TABLE  11.9 

SIMPLIFIED  NOZZLE  SIZING  PROCEDURE 


1.  NOZZLE  LENGTH  —  AS  LONG  AS  POSSIBLE  WITHIN  GEOMETRIC  CONSTRAINTS 
2  INLET  DIAMETER 


/ T - 

4  /ntj  \ 

P  g,in 

1  4»  rtf 

(p  v  ) 

'^g  g' in 

pi 

where  k s  Vg/Vj= constant 


(Eg  n  .12) 


3  THROAT  LOCATION  —  AT  50%  OF  NOZZLE  LENGTH 

THROAT  PRESSURE  IS  HALF  WAY  BETWEEN  INLET  AND  EXIT  PRESSURES 


4.  THROAT  DIAMETER 


■  m 


I  + ' 


;  c 


i 


thr  K  pi 

1  + 


g,in 


rP 


9, thr 


V- 


.  2RT,,m 

LP"  +1 

f  ^  P,hr^l 

9,«n  |+rK 

f  pthr  1 

(Eg  n.  13) 


,  .  /vg,in  w  pin  /* ~  6  in ~ \ 

thr=  InV  Vgi1hr  PthrX(l-ethr  ) 


(Eg  n.  15) 


5  EXIT  DIAMETER 


'exit  K, 


I  + 


r P  9, exit 


^.exit  "  V  vg,in  +  , 


2RTq,th  .  W.  Pex>t\ 

l  +  r*  in  Pexit  V-«in/\  "  Pin/ 


(Eg  n  .17) 


^fcxit  “  ^tn 


'  V 

vg,m 

P 

w  V 

vg,exit 

x  p  x 

Kexit  ' 

(Eg  II  18) 
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EFFECT  OF  EXIT  PRESSURE  ON  NOZZLE  PERFORMANCE  PREDICTS 
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TWO-PHASE 
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FIG  11.21 


COMPARISON  OF  MODEL  WITH  DROPLET  HEAT  TRANSFER  TO  EXPERIMENTAL 
WATER-NITROGEN  RESULTS  AT  NOZZLE  INLET  PRESSURE  OF  150  psia 
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EXIT  VELOCITY,  fl/sec 
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FIG  11.22 


COMPARISON  OF  MODEL  WITH  DROPLET  HEAT  TRANSFER  TO  EXPERIMENTAL 
WATER-NITROGEN  RESULTS  AT  NOZZLE  INLET  PRESSURE  OF  190  psia 
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EXIT  VELOCITY,  ll/sec 
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FIG  11.23 


COMPARISON  OF  MODEL  WITH  DROPLET  HEAT  TRANSFER  TO  EXPERIMENTAL 
WATER-NITROGEN  RESULTS  AT  NOZZLE  INLET  PRESSURE  OF  230  psia 


LOADING  RATIO.  mg/mg 


DIMENSIONLESS  NOZZLE  DIAMETER.  d/dm  NOZZLE  EFFICIENCY. 


9-4  FIG  II  25 

COMPARISON  OF  SIMPLIFIED  SIZING  PROCEDURE  TO  DROPLET  HEAT 

TRANSFER  MODEL 
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FIG  11.28 


SCHEMATIC  OF  WORK  EXTRACTION  SYSTEM  WITH  TWO-PHASE  NOZZLE 


heat  input 


heat  REJECTION 
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